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Abstract: If R-parity conserving supersymmetry exists below the TeV-scale, new par 
tides will be produced and decay in cascades at the LHC. The lightest super symmetric 
particle will escape the detectors, thereby complicating the full reconstruction of the decay 
chains. In this paper we expand on existing methods for determining the masses of the 
particles in the cascade from endpoints of kinematical distributions. We perform scans 
in the mSUGRA parameter space to delimit the region where this method is applicable. 
From the examination of theoretical distributions for a wide selection of mass scenarios it 
is found that caution must be exerted when equating the theoretical endpoints with the 
experimentally obtainable ones. We provide analytic formulae for the masses in terms of 
the endpoints most readily available. Complications due to the composite nature of the 
endpoint expressions are discussed in relation to the detailed analysis of two points on 
the SPS la line. Finally we demonstrate how a Linear Collider measurement can improve 
dramatically on the precision of the masses obtained. 
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1. Introduction 



The Standard Model (SM) of particle physics has been remarkably successful in describing 
the physics probed by modern day particle accelerators. No deviation from the SM has 
thus far been confirmed by experiment and only the Higgs mechanism, the SM's instrument 
for the breaking of the electroweak symmetry, remains to be discovered. Nevertheless, 
the SM suffers from considerable theoretical difficulties, not least of which is the hierarchy 
problem the extreme sensitivity of the electroweak scale to new physics. Such difficulties 
imply that the SM is only an effective low-energy theory (albeit a highly successful one) 
applicable only up to a few hundred GeV or so, and will need to be extended in order to 
describe physics at higher scales. 

One extension which has attracted a lot of attention is supersymmetry [||, [|, ^j, [5||. 
Supersymmetry not only solves the hierarchy problem but has many other attractive fea- 
tures: it is the only non-trivial extension to the Poincare symmetry of space-time Q; it 
is essential to the formulation of superstring theories 0; it provides a low-energy theory 
which is more amenable to the unification of the fundamental forces into a Grand Unified 
Theory (GUT) at some high energy scale ||; it provides a natural mechanism for generat- 
ing the Higgs potential which breaks the electroweak symmetry j|, [l0|, |ll], [12], 13, 14 1; and 
it supplies a good candidate for cold dark matter fl5[| . Furthermore, if it is to be relevant in 
solving the hierarchy problem it must exhibit experimental consequences at the TeV-scale, 
and therefore can be tested by experiment at the Large Hadron Collider (LHC). For an 
overview of supersymmetry searches at LEP, the Tevatron and HERA, see Ref. [fig] . 

If supersymmetric particles are produced at the LHC, thus confirming supersymmetry, 
it will become important to identify them and accurately measure their masses. This will 
be essential for identifying the low energy model and hopefully distinguishing the Minimal 
Supersymmetric Standard Model (MSSM) from other non-minimal extensions. Further- 
more, since no supersymmetric particles have so far been discovered, supersymmetry must 
be broken by some as yet unknown mechanism. Only an accurate determination of the su- 
persymmetric particle masses and couplings will allow us to determine the low energy soft 
supersymmetry breaking parameters. It is hoped that extrapolation of these masses and 
couplings to high energies using the renormalisation group equations will provide an insight 
into the mechanism of supersymmetry breaking and, more generally, physics at the GUT 
scale [17]. Since errors in the mass measurements will be magnified by the renormalisation 
group running it is absolutely essential that these masses be determined as accurately as 
possible. 

Here we will discuss supersymmetric mass measurements with reference to one particu- 
lar model of supersymmetry breaking, minimal super-gravity (mSUGRA) [18, 1^, ^ [l^, 14]. 
In this model, the supersymmetry is broken by the interaction of new particles at high en- 
ergy which are only linked to the usual particles by gravitational interactions; this new 
sector of physics is often referred to as the hidden sector. These gravitational interactions 
transmit the supersymmetry breaking from the hidden sector to our own sector, produc- 
ing TeV scale effective soft supersymmetry breaking terms in the GUT scale Lagrangian, 
quantified by parameters which run logarithmically down to the probed TeV scale. At the 



- 2 - 



GUT scale, the scalar supersymmetric particles are assumed to have a common mass, mo, 
while the gauginos have a common mass ra^j. The trilinear couplings are also taken to be 
universal at the GUT scale and denoted Aq. 

Mass measurements in the MSSM are complicated by R-parity conservation, which 
is introduced to prevent unphysical proton decay. R-parity conservation requires that 
supersymmetric particles are produced in pairs and causes the lightest supersymmetric 
particle (LSP) to be stable. Consequently the LSP is inevitably the end product of every 
supersymmetric decay and, if electrically neutral, will escape the detector leaving no track 
or energy deposit. While this provides a very distinctive missing energy signature, it makes 
it very difficult to measure masses at the LHC since one cannot fully reconstruct decays. 

Instead, mass measurements rely on continuous mass distributions of decay products 
which attain extrema for certain configurations of the particle momenta that are unam- 
biguously determined by the masses of initial, intermediate and final particles involved. 
These relations may often be inverted to give the masses of unstable particles. This is 
analogous to the way a bound on the neutrino mass can be obtained from the end-point of 
the beta-decay spectrum of 3 H [19], but is usually more complex, since a long decay chain 
is often involved. 

In this study we will consider supersymmetric mass measurements made by examining 
the mass distribution endpoints or 'edges' of the long decay chain 1 q — ► X2Q ~^ Ml ~^ Xi^Q 
in the ATLAS detector |p0| ]. In particular, we will focus on the Snowmass mSUGRA 
benchmark line SPS la [21], but will also include other mSUGRA parameters in a general 
discussion. In addition to the usual SPS la point, which we will denote SPS la (a), we will 
also consider another point on the SPS la line, denoted SPS la (/?), which has a reduced 
branching ratio for the decay. This will provide a counterpoint to the study of SPS la (a) 
where the branching ratio is rather high. This study differs from previous reports p^, |23], 



24, ^5|, by (i) discussing theoretical distributions which arise for different mass scenarios, 
(ii) providing inversion formulas, (iii) discussing ambiguities and complications related to 
the composite nature of the endpoint expressions. Furthermore, we provide an overview of 
the mSUGRA parameter space, and consider a new point on the SPS la line. Finally, we 
discuss the effects of including Linear Collider data in the analysis. 

In Sect. I we will discuss general mSUGRA scenarios, paying attention to supersym- 
metric decay branching ratios to gain some understanding of how generally applicable 
these kinematic endpoint measurements are. We will identify mass hierarchies for which 
these mass measurements are possible and show that they occur over a large portion of 
the mSUGRA parameter space. In Sect. || we will outline the properties and mass spec- 
tra of the Snowmass benchmark line SPS la and the associated points SPS la (a) and 
SPS la (/3). After defining the decay chain under investigation, we go on in Sect. I to 
discuss the theoretical framework of cascade endpoint measurements, and present analytic 
expressions for masses in terms of these endpoints. In Sect. || the experimental situation 
at the ATLAS detector will be studied, including event generation and reconstruction, the 
removal of backgrounds, and results for the measurement of the kinematic endpoints will 

1 Throughout the text the notions lepton/slepton as well as I /I will refer to the particles of the first and 
second generation. The third generation particles will be called tau/stau and denoted r/f. 



-3- 



be presented. The extraction of masses from the kinematic endpoints will be described in 
Sect. ||. Finally, the remarkable improvement of the accuracy of the mass measurements 
obtained by using inputs from an e + e~ collider [27, 28 1 will be studied in Sect. [7|, before 
drawing our conclusions in Sect. @. 



2. Cascade decays in mSUGRA scenarios 

In this paper we will be examining the decay chain q —* X2I ~^ ^1 ~^ Xi^Q i n the Snowmass 
scenario SPS la for the purpose of the extraction of the supersymmetric particle masses. 
However, it is extremely unlikely that SPS la is exactly the parameter choice of reality, 
and there would be little point to the study if our methods were only applicable at SPS la, 
or a small region around it. Therefore, in this section we will take a more general look 
at mSUGRA scenarios to determine whether or not these methods may be used more 
generally, over a wide parameter range. 

At the LHC the main supersymmetric production will be sparticle pairs gg, gq and 
qq, as long as these are not too heavy. Each sparticle immediately decays into a lighter 
sparticle, which in turn decays further, until, at the end of the chain, an LSP will be 
produced. Since there are two parent supersymmetric particles, each event will typically 
have two such chains, complicating their reconstruction. 

What route is taken from the initial gluino or squark down to the LSP depends on 
which decay channels are open, as well as their branching fractions. Therefore there are 
two criteria which must be met in order to use the cascade kinematic endpoint methods: 
firstly the sparticle mass hierarchy must be such that the analysed decay chain is allowed; 
and secondly, the cross-section for the entire decay chain must be large enough to allow 
analysis. 

2.1 The mSUGRA mass hierarchy 

In a general MSSM model there are few constraints on the sparticle masses, so little can 
be assumed about their relative mass hierarchy. However, if universal boundary conditions 
are imposed on the mass parameters at the GUT scale, some mass orderings at the TeV 
scale are natural, and some are even necessary. 

For example, in mSUGRA scenarios squarks and sleptons have a common mass at the 
GUT scale, mo, but when the masses are evolved down to the TeV scale, QCD interactions 
affect only the evolution of the squarks and not the colourless sleptons. Consequently, 
squarks are always heavier than sleptons at LHC energies. For the same reasons, although 
they may start off with different masses at the GUT scale, running induced by QCD 
interactions is usually enough to result in the squarks of the first two generations being 
heavier than the neutralinos and charginos at the TeV scale. (Large mixing can however 
bring t±, and to a lesser extent bi, quite low in mass.) 

At the opposite end of the spectrum, there are a number of possible candidates for 
the LSP: the lightest neutralino, Xii the lightest chargino Xi > or the lightest slepton (in 
mSUGRA models, the gravitino is usually rather heavy). Which of these is the LSP 
depends on the relative sizes of mo, m^j and the derived low-energy Higgs-higgsino mass 
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parameter fi. However, the assumption of gauge unification at the GUT scale explicit in 
mSUGRA models leads to the relation 

5 

Mi w - tan 2 6 W M 2 (2.1) 
3 

between the U(l) and SU(2) gaugino masses, Mi and M 2 respectively. As a result, Mi 
tends to be rather low, significantly lower than m\ji- Furthermore, the derived quantity 
fx is often required to be much larger than Mi in order to give the correct electroweak 
symmetry breaking (at the SPS la (a) reference point \x = 357.4 GeV). For the majority 
of parameter choices this implies that the LSP will be x?> with t\ being the LSP only if 
rriQ <C mi/2, and x^ only for a small region where — ► 0. The left-handed sneutrino, 
by virtue of its SU(2) interactions, is usually heavier than fx, and is anyway ruled out by 
direct searches [pfifl . It is indeed fortunate that Xi is the LSP for most of the parameter 
space since it is clear that only an electrically neutral LSP can play the role of the dark 
matter constituent which is believed to fill the universe. Finally, the gaugino mass relation, 



Eq. (2.1), implies that the LSP is usually bino-like. 

The first requirement for the decay chain q — > x^Q — * HQ ~^ Xi^Q is that the gluino 
should be comparable to or heavier than the squark initiating the decay chain. If the 
gluino is sufficiently light, then the squark will almost always choose to decay via its strong 
interaction q — > gq rather than by the electroweak decay q — > x^Q- Of course, one does not 
need all of the squarks to be lighter than the gluino; as long as one squark, for example b\, is 
lighter than the gluino, useful information can potentially be obtained from its subsequent 
decay chain. The second important characteristic is that x% should be heavier than I, 
thereby allowing the lower part of the chain to proceed, x§ ~~ * ^ ~~ * Xi^- Otherwise x\ 
will decay to Xi^ or XiK or to Xiff via a three-body decay, and the useful kinematic 
endpoints are lost. 

In order to understand where in the mSUGRA parameter space these hierarchy require- 
ments are realised, we have performed a scan over the my 2 -mo plane for four different 
choices of Aq and tan/3 (with fj, > 0), and identified the different hierarchy regions with 
different colours in Fig. |l|. The renormalisation group running of the parameters from the 
GUT scale to the TeV scale has been done using version 7.58 of the program ISA JET p0|j . 
which is inherent to the definition of the 'Snowmass Points and Slopes' (see Sect. ^). 

The upper left plot shows the mj^-mo plane with Aq = —uiq and tan/3 = 10 and 
includes the SPS la line and points (labeled (a) and (/3)). The upper right plot has Aq = 
and tan j3 = 30 and contains the benchmark point SPS lb. The lower left plot also has 
Aq = but tan /? = 10 and contains the SPS 3 benchmark line and point. Finally the lower 
right plot has Aq = —1000 GeV and tan/3 = 5 and contains SPS 5. 

The different hierarchies themselves are combinations of the hierarchy between the 
gluino and the squarks important to the upper part of the decay chain, and that of x® 
and the sleptons relevant to the later decays. Since rrt[ < mj L for any set of mSUGRA 
parameters, we here use Ir. The seven numbered regions are defined by: 

(i) g > max(di ,UL,h, h ) and xk > max([ R , fi ) 
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m, 3 [GeV] m, 3 [GeV] 



Figure 1: Classification of different hierarchies, labeled (i)-(vii), for four combinations of tan/3 and 
Ao, such that the four panels contain respectively the SPS la line, the SPS lb point, the SPS 3 line, 
and the SPS 5 point. The regions marked 'TF' are theoretically forbidden. (See text for details.) 



(ii) 


g > max(d£,tt£,6i,£i) 


and 


Ir > xl > h 


(iii) 


g > max(dL,ui J ,bi,ti) 


and 


min([ R ,fi) > xl 


(iv) 


d L > g> max(u L ,6i) 


and 


min^jfi) > xl 


(v) 


mm(d L ,u L ) > g >bi 


and 


min([ R ,fi) > xl 


(vi) 


mxa(dL,UL,h) > g > i\ 


and 


min(/ R ,fi) > xl 


Vii) 


mm.(d L ,u L ,bx,tx) > g 


and 


min([ R ,fi) > xl 



where for fermions a particle's symbol represents its mass, while for scalars a particle's sym- 
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bol represents the sum of the masses of the scalar and its SM partner. Also shown (mauve) 
is a region where the LSP is charged and therefore ruled out, as well as a theoretically 
forbidden (TF) region (gray) for low mi/2- 

It is interesting to note that there are no regions where a squark is heavier than the 
gluino and is heavier than one of the sleptons. This is simply because the gluino and 
[the gaugino part of] the neutralino have a common mass, m\/2i at the high scale, so if the 
gluino is light, the neutralinos will also be light. 

However, for more general non-mSUGRA unification scenarios one could still expect 
hierarchies of the type > rrig > m^o > rnj > m^o to be realised. It would then be 
important to be able to distinguish one hierarchy from the other; this should be possible 
using the kinematic endpoints, number of 6-quarks in the final state, etc. Also it is possible 

to distinguish xjj — ► M — ► X?^ ( m x? 2 > m l) fr° m X2 — * Xi^ ( m i > m x%)- ^ ne ^ rs ^ nas 
well-known triangular shape of mu while the second has a typical 3-body shape. All in all 
it should therefore be possible to distinguish the various hierarchies ( |2.2D before continuing 
to determine the masses themselves. 

Region (i) is the only one that has a 'useful' squark decay together with a decay of X2 
to a slepton, and is shown in light and dark green in Fig. [l]. We see that there is therefore a 
large region where the mass hierarchy is compatible with the methods presented here, and 
even though we will only perform the analysis for the points (a) and (/?) on the SPS la 
line, one would expect these methods to be widely applicable. 

However, these plots of the mass hierarchies really only show the regions in which the 
masses are such that the decay chain may occur. If the full decay chain is to be useful, 
it must have a sufficiently large branching ratio to be seen above the many backgrounds. 
We will therefore go on to examine the branching ratios of the pertinent decays over (a 
restricted range of) the mSUGRA parameter space. As a first taste, we have highlighted 
in a brighter green (and denoted (i) a ) the part of region (i) corresponding to where the 
overall branching ratio for the decay chain q — > X2I ~^ Mq ~^ Xi^Q exceeds a tenth of that 
at the SPS la (a) reference point. Although the decay chain is available over a rather 
large region of the parameter space, using this decay for large values of mi/2 and mo will 
be extremely challenging due to the small branching ratio. 

2.2 Gluino and squark decays: the upper part of the chain 

The decay branching ratios of the gluino are shown in Fig. ^ over the mx^—mo plane for 
two different scenarios. The representation is such that the branching ratio of a given decay 
channel in a small neighbourhood of the my 2 -mo plane is equal to the fraction which the 
corresponding colour occupies in that neighbourhood. Since the gluino only feels the strong 
force, it has to decay into a quark and a squark. If no squark is light enough, a three-body 
decay through an off-shell squark will take place; this is what happens in the green/white 
region of Fig. ^ at small Tn\i2- 

For the rest of the m^ 2 -wio plane the gluino decays fairly democratically into the 
accessible squarks. In considerable parts of the plane only one two-body decay is open, b\b 
(red) or tit (yellow), in which case the allowed decay takes close to the full decay width. 
Although one can in principle obtain information about the gluino mass by analysing its 
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Figure 2: Decay channels of g with = — mo, tan/3 = 10 (left) and Ao = 0, tan/3 = 30 (right). 
In the left panel the SPS la line is shown together with the two points (a) and (f3) marked with 
triangles. In the right panel the triangle marks the SPS lb point. The branching ratio of a given 
decay channel in a small neighbourhood of the m 1( / 2 -mo plane is equal to the fraction which the 
corresponding colour occupies in that neighbourhood. The region where Xi is n °t the LSP is 
denoted 'Charged LSP' and is discarded. Some regions are also forbidden theoretically, in that e.g. 
it is not possible to obtain electroweak symmetry breaking (labeled 'TF'). 



decay chain, we will only consider here the decay chain starting from a parent squark, and 
leave the gluino case for a separate publication [33|. 

As already intimated, squarks may decay by the strong force into a quark and a gluino 
(if the gluino is lighter), or decay by weak interactions into a quark and a chargino or 
neutralino, or via a loop into a gluon and a lighter squark. If kinematically allowed, the 
strong interaction takes a large fraction of the branching ratio, but since the (lighter) 
charginos and neutralinos are typically much lighter than the gluino, there will always be 
some neutralino production. 

Within mSUGRA models, the squarks di and ul are very close in mass and behaviour. 
Furthermore, the second generation squarks, §l and cl, are almost identical copies of the 
former two. It is therefore useful to have the common notation, q~L, for these four squarks. 
In a similar manner <jr is used for d,R, ur and their second generation copies. The right- 
handed squarks differ from the left-handed ones in that they do not feel weak interactions, 
which again makes their decay pattern different. In Fig. ^ the decay branching ratios of 
ul and ur are shown in the my 2 -mo plane for two different scenarios. 

For rrii/2 *C mo, when the gluino mass is smaller than the squark mass, both and 
even more so q~R have strong-interaction decays. In the rest of the m^-rao plane, when 
the strong decay is forbidden or suppressed by phase space, their decay patterns are very 
different: while q~R decays directly into the LSP, qi, prefers x° an d xf- 

For both low and high tan (3, the x± is predominantly bino, with only a tiny admixture 
of wino and higgsino, while x% an d xt are m ainly wino. For quite low mass parameters, 
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m 1iB [GeV] m 1i5 [GeV] 



Figure 3: Decay channels of uj, (top) and ur (bottom) with Aq — —mo, tan/3 = 10 (left) and 
Aq = 0, tan/3 = 30 (right). (See the caption of Fig. ||.) Note that the branching ratios for all first 
and second-generation squarks are assumed to be essentially the same. 

m i/2 ^5 100 GeV, they become more mixed. Since qz, generally has a much larger SU(2) 
coupling than U(l) coupling, decays to X\ an d X2 wm be preferred unless the difference in 
phase space makes the decay to the lighter Xi competitive. In contrast, since the (Jr has 
no SU(2) interaction it will decay predominantly to the bino X\i except at quite low mass 
parameters where the neutralinos change character. 

The third generation squarks differ from the others in two aspects. First, the mass 
eigenstates can have more even admixtures of both handedness components. For b this is 
the case for low mass parameters, mo, mi/2 ^ 200 GeV, where the branching ratios into Xi 
and X2 are °f comparable size. At higher masses 61 ~ bi and 62 ~ t>R, giving a b\ which 
prefers to go to X2 rather than Xi- Second, due to large splitting, the third generation 
squarks can decay into other third generation squarks together with a weak gauge boson. 
The drastic change observed in figure || as rno,mx/2 become less than ~ 200 GeV, is due 
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Figure 4: Decay channels of b\ (upper panels) and 62 (lower panels) with Aq — —too, tan/3 = 10 
(left) and A = 0, tan/3 = 30 (right). (See the caption of Fig. f§) 



both to the more mixed mass eigenstates for lower masses, and to the closing of certain 
channels involving t or a heavy gauge boson. 

While b% has a large branching ratio into X2 throughout the entire plane, 62 produces 
X2 a t a much smaller rate, except for small mass parameters. 

To summarize, the squark decays that are 'useful' for kinematic endpoint analyses, 
are those of left-handed first and second-generation squarks, as well as those of b\ and to 
a lesser extent 62- These occur in the entire mi/2-mo plane, except for extreme values 
mp2 <C mo, and for both low and high tan/3 values. For quite low mass parameters also 
qjl contributes. 

2.3 Neutralino and slepton decays: the lower part of the chain 

The reasons why X2 often plays an important role in the reconstruction of SUSY events are 
many. Kinematically situated midway between the initially produced gluinos/squarks and 
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the LSP, it is abundantly decayed into, as we have seen. What makes it so useful, usually 
more so than the x^, which is produced in similar ways and amounts, is the fact that its 
decay products, in addition to easily setting off the trigger, also reconstruct well. 

In Fig. [| the main decay channels of X2 are shown for two values of tan j3. A two-body 
decay is preferred over a three-body decay, and the coupling to Z, f or v is usually stronger 
than the coupling to Xi- 

For mo > 0.5m l / 2 all of I, f and v are heavier than \\i so oru y the decay into the LSP 
is possible. In yellow, to the very left, no two-body channel is open, and X2 undergoes a 
three-body decay, proceeding through an off-shell squark or stau/slepton, or involving an 
off-shell Z, W or h. As m^j increases, also m^o — m^o increases and more decay channels 
become available. First the Z channel opens and takes the full decay width, then the h 
channel opens to dominate. The mass difference between x? and X2 i s mostly independent 
of mo, which is why, to a good approximation, the yellow, red and green regions are stacked 
horizontally. 

In the blue regions decays into l/f/v are kinematically allowed. Following a clockwise 
movement, the scalar masses are reduced relative to x<]- The right-handed scalars are 
lighter and become available first (dark blue region). In the light blue region the left- 
handed scalars have become available and, despite less phase space, take most of the width 
due to their SU(2) coupling. The black part of the blue regions shows the decay into Ir 
(dark blue region) and II (light blue region). These are the decays of interest to us. At 
low tan/3 (left panel) the slepton channels can be used in most of the blue regions. At high 
tan/3 the situation is less optimistic. The fi channel totally dominates the Ir channel, and 
only in a small region is the II channel open. 

Decay products which involve tau particles are more difficult to use since their recon- 
struction is always incomplete due to undetected neutrinos. However, in some parts of the 
parameter space, especially at high tan (3, these channels take the full decay width, so one 
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must be prepared to use them. In any case, whatever the MSSM parameters should turn 
out to be, in order to measure the stau mass itself, tau particles must be reconstructed. 

In the non-blue regions of the plane, information on the sparticle masses can still be 
retrieved, although to a lesser extent. From an experimental point of view final states 
involving two leptons are preferable. In the yellow region this fraction is marked in black. 
In the red region Z decays leptonically in 7% of the cases. In the green region, no leptonic 
decay is available. Here the bb final state of the Higgs can be used. This channel may even 
serve as a discovery channel for the Higgs boson. 

Returning once more to our chosen decay chain, q — > x®Q ~^ MQ ~~ * Xi^9> ^ is clear 
from Figs. that the initiating squark must be or b. Furthermore Fig. |H| shows that 
for low tan/3 the sleptonic decay of X2 ^ s open in a large fraction of the my 2 -mo plane. 

2.4 Other constraints 

Much of these parameter planes considered in Figs, [ij — |sj are actually excluded or disfavoured 
by observations. The left part (small rni/2) is typically excluded by the lower bound on 
the Higgs mass, and the region of 'large' mo + m\j2 is excluded by the WMAP data 
p9| , |40p ) since a too high contribution to the Cold Dark Matter (relic LSP) density is 
produced. Such bounds have been explored in considerable detail for the so-called 'Post- 



LEP' SUSY benchmark points |35|]. For the SPS points, see Ref. [36|. Also, there are 
constraints from the non-emergence of unphysical vacua during the renormalisation-group 
running from the high scale |$7|]. These typically rule out a sector at low mo and high 
mi/2 that may extend beyond that excluded by charged LSP. We shall here ignore such 
additional constraints, since they depend somewhat on the assumptions which are adopted. 
In particular, tunneling into them may take longer than the age of the Universe. 



3. Summary of SPS la 

We have shown that the squark and gluino initiated cascade decays seen at the benchmark 
point and slope SPS la are not atypical of a large portion of the parameter space. In this 
section we will go on to explore SPS la in more detail, and introduce our second SPS la 
point on the line. 

3.1 The SPS la line and points 

The SPS la benchmark line is defined as the masses and couplings of supersymmetric 
particles at the TeV scale as evolved from the GUT scale mSUGRA inputs 

m = -Aq = OAmi/2, 

tan/3 = 10, H>0, (3.1) 



by version 7.58 of the program ISA JET [3C]. Elsewhere in this report when the mSUGRA 
GUT scale parameters are referred to it is to be understood that the low energy parameters 
are obtained in this way. 
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Figure 6: Masses (left) and relative widths (right) of relevant sparticles as m\/2, mo and Aq are 
varied along the SPS la slope, defined by Eq. (EO). The vertical dotted lines represent SPS la 
points (a) and (/?). 

In addition, we define two points (a) and (/?) on the SPS la line according to: 
(a) : m = 100 GeV, m l/2 = 250 GeV, 

(/?) : m = 160 GeV, m l/2 = 400 GeV. (3.2) 



The first point (a) is the 'basic' SPS la point of Ref. [21] and studied in Ref. [27, 28|, while 
the second (/3) is a new, less optimistic scenario with a reduced cross-section for the decay 
chain. 

The masses of particles relevant for our analysis are shown in Fig. |6| (left), moving 
along the SPS la line by varying m^j. The values of the masses at points (a) and (/?) 
can be seen from the vertical dotted lines. As expected, all the masses except the lightest 
Higgs boson mass increase linearly with mi/2- Neither the heaviest neutralino mass nor the 
chargino masses are shown; to a good approximation, m~± ~ m^o ~ m^o and m~± — m^o. 
Similarly, the masses of and A are not shown, but m^± ~ uia — ran- 

For the points (a) and (/3), these masses are further detailed in Table [l], with the 
masses of the particles in our chosen decay chain displayed in bold. 

The relative widths (width divided by the mass) of the decaying sparticles are shown 
in Fig. H (right), and are everywhere less than 2% of the mass. The wiggles in some of 
these curves (as well as in some of the branching ratio curves below) are due to limited 
precision in ISA JET |38|j . As will be discussed, these widths contribute to a blurring of the 
kinematical endpoints, and will thus be reflected in the mass determination. 

3.2 Sparticle production 

The cross-sections for producing supersymmetric particles at the LHC are for moderate 
values of m i / 2 rather high. This can be seen in Fig. |7| which shows the dominating sparticle 



- 13 - 



Point 


Q 




d>R 


ut 




b 2 


h 


h 


ti 


(a) 


595.2 


543.0 


520.1 


537.2 


520.5 


524.6 


491.9 


574.6 


379.1 


(P) 


915.5 


830.1 


799.5 


826.3 


797.3 


800.2 


759.4 


823.8 


610.4 




e T 


^ it 


To 


Tl 




Vt 

T L 






A 


(a) 


202.1 


143.0 


206.0 


133.4 


185.1 


185.1 




401.8 


393.6 


(P) 


315.6 


221.9 


317.3 


213.4 


304.1 


304.1 




613.9 


608.3 




xl 


xl 


xl 


X? 


X.2 


xt 




H 


h 


(a) 
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114.0 
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538.4 


299.1 


161.0 


553.3 


299.0 




608.9 


117.9 



Table 1: Masses [GeV] for the considered SPS la points (a) and (/?) of Eq. (3.2) 




pair production cross-sections, as m^ 2 is varied along the SPS la line. Notice that these 
cross-sections fall very rapidly as m 1( / 2 is increased, which will cause repercussions in the 
analysis of SPS la (/3). 

The cross-sections for gluino-gluino, gluino-squark and squark-squark pair produc- 
tions are detailed in Table § for the two chosen analysis points, together with the SUSY 
total rate. Of course since other supersymmetric particle pairs may contribute to the total 
SUSY rate it is not simply a sum of the other numbers in the table. 

These supersymmetric particle pairs are predominantly produced by QCD interactions 
of quarks and gluons in the colliding protons. For gluino pairs this is mainly due to gg — ► gg 
via t-channel gluino exchange and s-channel gluons, and at a much smaller rate qq — > gg via 
s-channel gluons. Squark pairs with the same handedness have the dominant production 
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process qq — * qL/R Q.L/R with a t-channel gluino exchange, but can also be produced via 
99 ~^ Ql/rQl/r with an s-channel gluon, or via qq — ► <Jl/r Ql/r with both t-channel 
gluino and s-channel gluons. For opposite handedness, qLQR production, the s-channel 
gluon exchange processes [and thus gg — > q^ qn at tree-level] are disallowed. Lastly, the gq 
final states are produced by the process gq — > gq mediated by t-channel squark or gluino 
exchanges and s-channel quarks. 





a(SUSY) 


°{gg) 




°{Mr) 








(a) 


49.3 


5.3 


11.4 


12.3 


3.5 


4.8 


4.1 


(A 


4.76 


0.29 


0.97 


1.06 


0.44 


0.61 


0.53 



Table 2: Selected supersymmetry cross-sections in pb. 



However, these particle pair cross-sections are not the production cross-sections that 
are relevant to our analysis of the decay chain, since it does not matter from where the 
parent squark originates. Therefore we should also be counting, for example, gluinos which 
decay into squarks as possible sources of the decay chain. In Table |3| we divide the sparticle 
productions rates into 'direct' and 'indirect' contributions, reflecting production rates from 
sparticles in the 'initial' supersymmetric state as opposed to their generation from the 
decay of a parent sparticle. Furthermore, if the 'initial' supersymmetric state contains 
two possible parents then the chance of generating the desired decay chain is doubled 2 . 
Therefore the pair-production rates with two possible parent particles are counted twice. 
Again, one cannot simply add the various contributions from Table ^ to obtain the 'direct' 
rates of Table ||, since they include contributions from squarks or gluinos produced in 
association with other supersymmetric particles. 
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0.25 


0.64 


1.44 




Total 


2.71 


3.21 


3.79 


0.50 


0.31 


0.87 


1.67 



Table 3: Selected sparticle production rates in pb. 



For low m-i/2 values Wln decay dominantly to a stau and a tau, see Fig. |8[ Here, the 
decay mode of interest is to a right-handed slepton, Ir, and a lepton, and is at the level of 
10%. For higher values of m\/2, the mass difference between Xi an d x% grows sufficiently 
to allow the decay \\ — ► y§h. At (/?), even though the Higgs channel takes a significant 
30% of the decay width, it was not used in this study. Due to the accuracy of lepton 
reconstruction compared to jet reconstruction, only in the case of very low lepton channel 
statistics can the Higgs channel improve on the results obtained with the slepton channel. 

2 The total branching ratio for the decay is sufficiently small that the chance of generating two of the 
desired decay chains is tiny. 
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Figure 9: The SPS la cascade decay chain. 

3.3 The cascade 

The cross-sections and branching ratios of our chosen decay chains 



Ql - 








(3.3) 


61 - 




- iptb - 


- Xllfltb 


(3.4) 


62 - 




- tpib - 


- Xilptb 


(3.5) 



are summarized in Fig. || for the two SPS la points. Since the left-handed up and down 
squarks, ul and d^, have very similar masses (at 537.2 GeV and 543.0 GeV respectively), 
for (a), they are in the above jointly referred to as gx,, and for this analysis will be grouped 
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together. For the fraction of decay chains which commence from a sbottom, b\ is responsible 
for 78% or so, leaving us rather insensitive to the contribution from 62- 

4. Mass distributions 

The longer a decay chain is, the more information it contains. To extract the masses of 
the supersymmetric particles in the decay we require at least as many kinematic endpoint 
measurements as unknown masses. In the lower part of the decay chain, where the second- 
lightest neutralino decays via x% ~^ W — * Xi^j there are three unknown masses: m^o, m i R 
and m^o. However, only two particle momenta are measured, those of the two leptons, 
from which only one mass distribution can be constructed, mj;. The system is highly 
underdetermined; one cannot extract the three masses, only a relation between them. 

When a squark is added to the head of the decay chain, q — > X2Q ~^ Wn<7 — * XiklnQ, 
three particles can be collected, and one can construct four mass distributions, mu, m q i n , 
m q i t and m q u, where following the notation of Refs. [El], 32 1, we denote the first emitted 



lepton Z n ('n' for 'near') and the second Zf ('f for 'far'). In principle this is just sufficient 
for extracting the four unknown masses: 1715, m^o, mr and m^o. However, in order to use 
the distributions m q i n and m q i { , we need to be able to distinguish Z n from Zf. Since this is 
usually not possible, two alternative distributions are defined, m^high) and m g ui ow ) |pl| , 
constructed by selecting for each event the largest and smallest values of m q i respectively. 

As will be detailed later in this section, the expressions for these kinematic endpoints 
are not always linearly independent, so these four endpoints are not always sufficient to 
determine the masses in the decay chain. In this circumstance one must look for other 
endpoint measurements. 

Correlations between different mass distributions can provide further measurements. 
For example, one may define the mass distribution mqn^z) identically to the m q u distri- 
bution but with the additional constraint 



mr7V2 < m u < m^. (4.1) 

This cut on mu translates directly into a cut on the angle 6 between the two leptons in the 
rest frame of Ir jO]. In terms of this angle, mjj is given by 



mil 



V(l -cosfl)/2 (4.2) 



so a constraint of the form (4.1) directly corresponds to 6 > |. The simplicity of this 
constraint allows one to find an analytic expression for the minimum of the mqii(Q>iL) 
distribution. 

In principle, other correlations between mass distributions could be used, but they are 
limited by the lack of analytic expressions for the associated extrema. It is no doubt possible 
to construct simple constraints for which analytic expressions for minima or maxima of mass 
distributions are possible, but this will not be investigated further in this study. 

If we were to also include a parent gluino at the head of the decay chain, g — > qq — > 
X2QQ ~^ IrIqQ ~^ Xi^QQi we would have an extra quark momentum at our disposal and 
could construct with it seven more ('primary') mass distributions. 
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4.1 Theory curves of invariant mass distributions 

In Fig. |l0| we show 'theory' versions of the five mass distributions discussed above for SPS la 
(a) and (/?), and three other mass scenarios. These distributions reflect the parton level 
only, where the quark and leptons are assumed to be perfectly reconstructed, and particle 
widths have been neglected, suppressing a mild smearing of the distributions. Leptons and 
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quarks are assumed massless; at worst, i.e. for an endpoint involving a 4.8 GeV 6-quark, this 
approximation gives a value for rn^ x which is wrong by 5.7 MeV at SPS la (a), which 
is negligible. Furthermore a common squark mass is used; in reality, the experimental 
distributions would be a sum of several similar distributions shifted typically by a few 
percent, depending on the differences between the squark masses. 

These theory curves have been generated from the phase space only, with no matrix 
elements inserted. For the two-body decay of scalar particles this is not a problem; the 
daughter particles are always emitted back-to-back in the rest frame of the parent and 
since a scalar particle provides no intrinsic direction the decay will be isotropic. The mu 
distribution, for example, is very well described by phase-space considerations alone since 
the scalar In removes any spin-correlations between the two leptons. For the decays of 
fermions this is not so clear; the parent's spin picks out a particular direction so the decay 
need not be isotropic. One would expect that spin-correlations must be fully taken into 
account by calculating the full matrix elements for the decay chain. In practice, while 
these spin-correlations are indeed very significant in individual processes, when one sums 
up final states containing leptons with positive and negative charge (which are produced in 
equal proportion) these spin correlations cancel out [ 43 1 . We have confirmed this assertion 
analytically. 



Most of the distributions of Fig. |10| show a strong dependence on the scenario. Only 
the mu distribution is independent of the masses involved, which is easily understood from 
( |4.2| ) together with the isotropic decay of Ir. The shape is easy to fit, and since only leptons 
are involved, the experimental resolution is high. The presence of the quark reduces the 
precision with which the other distributions may be measured. For m q u there is a clear 
dependence on the masses involved, the main feature being the length of the straight 
section, ranging from infinitesimal to the entire distribution. However, the maximum edge 
itself is quite well described by a straight line, at least for the latter part. Next, the 
constrained m q u distribution differs from m q u for smaller invariant masses. Its rise from 
threshold is not very well described by a straight line making the measurement of the 
minimum value rather imprecise. Both m q i distributions are 'composite' distributions, 
based on entries from m q i n and m q i f . This double nature is readily apparent for m g z(high)> 
which has a wide variety of shapes. Particularly dangerous are mass scenarios like {(5) or 
(ii), where the 'foot' which forms the last part of the edge can be hidden by backgrounds, 
giving a false maximum. In Sections ||-|6| these considerations will come up in relation to 
(/3). For m g ;(i ow ) a similar danger arises in cases like (i) where what appears to be an 
approximately linear descent to zero suddenly turns into a vertical drop towards the end. 

These theory curves must serve as guidance for our determination of the endpoints. 
As a first stage they should be the basis for the choice of functions with which the edges 
of the different distributions are fitted. As a second stage the more ambitious goal may be 
to fit not only the endpoint, but the entire distribution. 

4.2 Formulae for kinematic endpoints 

The invariant masses of various subsets of particles can be determined from kinematical 
endpoints and thresholds, as discussed in [MJ. We confirm the results relevant for our 
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analysis: 
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where 'low' and 'high' on the left-hand side in Eq. (4.5) refer to minimising and maximising 
with respect to the choice of lepton. Furthermore 'min' in Eq. (|4.9|) refers to the threshold 
in the subset of the m q u distribution for which the angle between the two lepton momenta 
(in the slepton rest frame) exceeds vr/2, corresponding to the mass range (|4.1| ). 

Notice that the different cases listed in Eq. (|4.4f) are distinguished by mass ratios of 
neighbouring particles in the hierarchy, m^/m^o, m^o/mj R and m^/m^o. Since each 
decay in the chain involves two massive particles and one massless one, the boosts from 
one rest frame to another are conveniently expressed in terms of such mass ratios. 
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4.3 Inversion formulae 

Once the endpoints of the distributions have been measured, the masses may be extracted. 
In principle, this can be done in two ways. If the number of endpoint measurements 
coincides with the number of unknown masses, one may analytically invert the expressions 
for the endpoints to give explicit formulae for the masses in terms of the endpoints. If more 
endpoints are known, the system becomes overconstrained and the measurements must be 
weighted according to their uncertainties. The analytic method cannot easily handle such 
a situation. Instead the more flexible approach of numerical fit must be used. 

However, the analytic inversion is more transparent than the numerical treatment, 
and reveals some interesting features of how the masses are related to the endpoints. It is 
also practical to use the analytical method in combination with the numerical method to 
provide initial values for the fits. Consequently we will discuss the inversion of the endpoint 
formulae in some detail in this subsection. 

The four principal endpoints ?7i™ ax , rn^ x , m ™^ w j and rn q ^^ is ^ are given by Eqs. (|4.3| )- 
( |4.5| ). While m™ ax is given by one unique expression, the others have different represen- 
tations for different mass spectra. To perform the inversion, each combination of the end- 
point expressions must be considered separately. There are four representations of rn^f* 
and three representations of (w^g*^ > m ^/(hi g h) ) ■ Each overall combination corresponds to 
a unique region in mass space (mq L ,rn^o,rnf R ,rn^o) and can be labeled by (i,j), where i 
and j denote the region of applicability of the m q u and m q i endpoints, respectively, as given 
in Eqs. Q-Q. 

However, not all of the 4x3 combinations are physical. In particular the regions (2,1), 
(2,2) and (3,3) are not possible. For (2,1) this can be seen by simultaneously trying to 



satisfy the mass constraints (2) of Eq. (|4.4| ) and (1) of Eq. (4.5). 

Additionally, for the regions (2,3), (3,1) and (3,2) the endpoint expressions are not 
linearly independent. In each of these three regions, one of the mass ratios is dominant, 
and a lepton (rather than a quark) is emitted from the vertex involving the dominant mass 
ratio. The maximum of m q u occurs when this lepton travels in one direction and the other 
lepton together with the quark travels in the opposite direction, all observed from the rest 
frame of any one of the sparticles in the decay. For massless quarks and leptons one can 
always write 

™%U = 2 PqPln + 2 PqPk + 2 PlnPh = m J(high) + m2 ql(low) + m l ( 4 - 10 ) 

In the three regions mentioned above, for maximum values this equation reduces to 

(m^n 2 = K^igh)) 2 + M") 2 (4.11) 

Since one of the endpoints can be expressed in terms of two others, the four endpoint mea- 
surements only transform into three independent conditions. Hence, in order to determine 
the masses an additional endpoint would be required for these three regions. The endpoint 
of m q ii^ g> 2L) is therefore particularly important in such cases. However, the complicated 
nature of Eq. (|4.9| ) leads to exceptionally cumbersome inversion formulae and will not be 
studied here. 
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For the remaining six regions the inversion of the four endpoint expressions Eqs. (|4.3| )- 
( |4.5| ) is possible. Coincidentally these are actually the regions most likely to be realised in 
SUSY scenarios with sparticle mass unifications at the GUT scale. Inversion formulae for 
these regions are detailed below. For simplicity of notation we will write: 
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rrig L = m^o + b (4.28) 
Region (4,2): 
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c 2 
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m ? - L = m^o + 6 (4.36) 

For a given set of endpoint measurements the inversion formulae of most regions will 
return unacceptable masses. In some regions the masses returned will be in contradiction 
with the presupposed hierarchy, rriq L > m^o > raj- > m^o, in some there will even be 
negative or imaginary masses. Such solutions must be discarded. If the masses returned 
by the inversion formulae of a region (i,j) do obey the hierarchy, they constitute a physical 
solution if they themselves belong to Solutions which do not satisfy this last constraint 
will be referred to as 'unphysical'. In principle it should be possible to construct for each 
region a number of conditions on the endpoint values which, beforehand, would tell if the 
inversion formulae would return a physical solution or not. Such conditions will however 
become quite complicated, and have not been sought in this study. Instead the more 
straightforward approach is taken: try the inversion formulae of each region, then discard 
any unacceptable solutions. 

Originating from realisation (4) of m^ ax , see Eq. ( [4.4] ) , where masses appear unsquared, 
the inversion formulae of region (4,3) in principle come in two versions, the difference being 
a sign within one of the mass expressions. In regions (4,1) and (4,2) it turns out that only 
one of the signs returns masses in the appropriate region. In Eqs. ( [4.27 ) and ( [4.29| ) the 
physical sign is therefore chosen. For (4,3) both versions must be considered. 

Once the required endpoints are measured and the inversion formulae used to determine 
the masses, we encounter a delicate feature of the entire method of obtaining masses from 
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endpoints. While the endpoints are given by single-valued functions of the masses, albeit 
with different expressions for different mass regions, the inverse is not true. A given set of 
endpoint values can in principle correspond to several sets of mass values. This is equally 
true for the numerical fit method, and has not received much attention previously (see, 
however, Ref. |32|]). This complication will be faced in Section ^. 



5. 'Data' generation and reconstruction 
5.1 Event generation 

The SPS points are defined by the low-energy MSSM parameters produced by ISAJET 7.58 
|30f|, given a set of high-energy input parameters. In our analysis PYTHIA 6.2 [[II] with 



CTEQ 5L 1 45] is used to generate the Monte Carlo sample. 3 To allow for this the low-energy 
parameters from ISAJET are fed into PYTHIA via the standard interface. PYTHIA in 
turn calculates the decay widths and cross-sections. Each event produced is passed through 
ATLFAST 2.60 @, a fast simulation of the ATLAS detector. In ATLFAST the output 
particles of PYTHIA are mapped onto a grid of 'calorimetric cells' of a given spacing 
in pseudorapidity rj and azimuth angle (j). Next, the cells are combined into clusters, 
after which particle identification takes place, including smearing of the four-momenta 
according to particle type. Jets are built by a cone algorithm with AR = 0.4, where 
AR = {Arf) 2 + (Acj)) 2 . Acceptance requirements are imposed: \i]\ < 2.5 for e/fi and 
\t]\ < 5 for jets as well as px > 5/6/10 GeV for e/^/jets. Leptons are marked as isolated if 
there is no other cluster within a distance AR = 0.4 of the lepton, and if additional energy 
inside a cone AR = 0.2 is below 10 GeV. While ATLFAST captures quite well the main 
features of the full simulation, some important effects are left out. Lepton identification 
efficiencies are not parametrized. A conventional 90% efficiency per lepton is therefore 
included by hand in the analysis. Also, the possibility of misidentifying a jet as a lepton is 
absent in the fast simulation, and has not been included in our analysis. The effect of pile- 
up on the jet energy resolution is accounted for in ATLFAST when run at high luminosity, 
as in this analysis, but pile- up events have not been simulated, and the underlying event is 
probably too 'slim'. However, as the selection criteria on jets and leptons are quite hard, 
we do not expect a more realistic detector simulation to change the results very much. 
Nevertheless, the numbers quoted at the end of this section should be validated with these 
effects included. 

The signature of a signal event is two opposite-sign same-flavour (SF) leptons, consid- 
erable missing pt from the escaping LSPs, and at least two hard jets, one from the signal 
chain, the other from the decay of the squark nearly always present in the other decay chain. 
The most important Standard Model process to have the same features as the signal, is tt 
production. Also W/Z together with jets, one of which is a 6-jet, can mimic the signal, 
and in combination with the underlying event, pile- up and detector effects, other processes 
will also now and then result in the given signature. Together with tt we therefore include 
QCD, Z/W+jet as well as ZZ/ZW /WW production. No if-factors have been used. 



3 The main parts of the analysis have been confirmed with HERWIG 
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The precuts (not tuned) used to isolate the chain are the following, 



• At least three jets, satisfying: pj, > 150, 100, 50 GeV 

> max(100 GeV,0.2M efr ) with M eS = E TMss + ELi^Ti 

• Two isolated opposite-sign same-flavour leptons (e or fx), satisfying p^P > 20, 10 GeV 

The QCD background is cut away by the requirement of two leptons and of considerable 
missing px- For the processes involving Z and W the requirement of high hadronic activity 
together with the missing px removes nearly all events. After the rather hard cuts listed 
above, the Standard Model background consists of approximately 95% ti. 

5.2 Different Flavour (DF) subtraction 

At this stage the main background events, in addition to ti, will come from other SUSY 
processes which have two opposite-sign same-flavour leptons. It is useful to distinguish 
between two types. The first, which will also be referred to as 'lepton-correlated', produces 
correlated leptons, e.g. leptonic decay of Z. In these processes the leptons always have 
the same flavour. The other type produces uncorrelated leptons which need not be of the 
same flavour. Typically the uncorrelated leptons are produced in different decay branches 
within the same event. Lepton universality implies that electrons and muons are produced 
in equal amounts (apart from negligible mass effects). This means that for events which 
produce uncorrelated leptons, there should be as many opposite-sign different-flavour (DF) 
lepton pairs produced as there are opposite-sign same-flavour lepton pairs, and their event 
characteristics should be the same. The same-flavour leptons are already part of the selected 
events. If one then subtracts the different-flavour events from the total same-flavour sample, 
this corresponds statistically to removing the non-signal same-flavour events which come 
from uncorrelated leptons. The residual of the subtraction is a larger statistical uncertainty 
in the new signal distribution. Clearly the different-flavour subtraction is a very effective 
tool which takes care of both ti and most of the SUSY background. Only processes with 
correlated leptons remain. 

In a more detailed study one would need to correct for the different acceptances of 
electrons and muons as functions of (px,rj,4>) m t ne different-flavour subtraction. Here we 
have assumed equal acceptance for electrons and muons, which even for a fast simulation 
is only approximately true. Also we have assumed that the efficiencies for reconstructing 
two same-flavour and two different-flavour leptons are the same. For close leptons this is 
not the case in the ATLFAST reconstruction. This will have some visible effects in the 
different-flavour-subtracted distributions, see below. 



Fig. 11 shows the different-flavour subtraction for the mu (top) and rn q u (bottom) 
distributions at points (a) (left) and (/?) (right). (The jet used in the m q u distribution 
is the one of the two pr-hardest jets which gives the smaller value of m q u.) All the plots 
have the same colour code. In red the same-flavour ('SF') distribution is shown by solid 
curve, the different-flavour ('DF') one by dashed curve. Their difference, the different- 
flavour-subtracted distribution ('SF-DF'), is shown in black with errors. The curve in blue 
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Figure 11: Different-flavour subtraction for mu and m q u at (a) (left) and ((3) (right). The 
solid/dashed red curves are the same-flavour ('SF')/different-flavour ('DF') distributions. In black, 
their difference, the different-flavour-subtracted distribution ('SF-DF'), is shown with error bars. 
The blue curve shows the part of the subtracted distribution which contains a signal chain ('SC'). 
The solid/dashed green curves ('SM') give the Standard Model part (completely dominated by tt) 
of the same-flavour/different-flavour distributions. They are statistically equal and will cancel each 
other. 



shows the part of the different-flavour-subtracted distribution which contains a signal chain 
('SC'). The reason why the blue distribution of m q u does not have the form of the theory 
distribution in Fig. is that the jet is only correctly selected in roughly half of the cases. In 
solid/dashed green the Standard Model contributions to the same-flavour /different-flavour 
distributions are shown ('SM'). They are statistically identical and will cancel each other 
through the different-flavour subtraction. 

If the samples contained only background events with uncorrelated leptons (and the 
different-flavour-subtraction procedure removed all of these) , the different-flavour-subtrac- 
ted distribution should fall exactly on top of the blue line. When this does not happen, 
it implies that the sample also contains background events with correlated leptons. The 
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Z-peaks in the ran distribution at both (a) and (/?) are obviously of this type. The Z's 
stem predominantly from the decay of heavier gauginos into lighter ones. At (a) ~ 80% 
of the peak comes from the heavy gauginos X3, x\ and % 2 ■ At (/?) X2 ^ s sufficiently heavy 
to decay into \\Z and is responsible for ~ 40% of the peak. The rest involves the heavy 
gauginos. 

For mu i$ 20 GeV a bump is visible in the different-flavour-subtracted distribution. 
This excess turns out to come mainly from lepton-uncorrelated events, predominantly of 
the type \\ — * which is abundant in the two scenarios, and where the taus decay 
leptonically. Since a pair of taus produce same-flavour and different-flavour leptons in 
identical amounts, the different-flavour subtraction should take care of this background. 
When it does not, it is because of an asymmetry in the reconstruction algorithm of ATL- 
FAST which accepts close same-flavour leptons at a higher rate than close different-flavour 
leptons. Since small AR between the leptons means small mu values, the bump appears 
at low mu. Such a reconstruction asymmetry may also have noticeable effects in other 
distributions, e.g. m q u, for which the maximum values in our two mass scenarios appear 
for parallel leptons. A dedicated study of the impact detector effects of this type may have 
on the endpoint determinations may be worthwhile. 

In the m q ii distributions there is an excess of events (compared with the blue line) 
at lower masses. Some of this excess comes from the lepton-uncorrelated events in the 
mu bump. Typically these land at low m q u values. The main contribution to the excess 
is however from xj^'s which decay sleptonically, but which originate from qji, t±, g or the 
heavier gauginos. 

At (/3) two effects complicate the endpoint measurements. One is a reduction of the 
SUSY cross-section by one order of magnitude relative to (a) . This allows for the Standard 



Model background to have a larger impact, see the green curves in Fig. 11. Since this 
background (practically only tF) consists of uncorrelated leptons, it is dealt with by the 
different-flavour subtraction. It leaves however an increased statistical uncertainty in the 
resulting distribution. The other effect at (/?) is a reduction in the signal branching ratio, 
mainly from a reduced BR(x2 — * IrI), see Fig. [| As a result also the SUSY background 
becomes larger relative to the signal. 

While the Standard Model background is practically negligible at (a), it is of similar 
size as the SUSY background at (/3). For the m q u distribution both background types are 
important only at low values, a fair distance away from the kinematical maximum value. 
The same is true for the m q i distributions. The endpoint measurements are therefore only 
minimally affected by the background. An exception is the mgu^y^ distribution, where a 
minimum is measured. 



5.3 Selection cuts 

A fair fraction of the events which pass the precuts, ~ 30% for (a) and ~ 20% for (/?), have 
more than two leptons. Each opposite-sign lepton pair which satisfies pt > 20, 10 GeV is 
used, either in the same-flavour or the different-flavour distribution. As usual, the different- 
flavour subtraction takes care of the incorrect combinations. 
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Figure 12: Quark selection cuts for m q u (left) and ra g m ow ) (right). Left: The two lower black 
curves in solid and dashed show the rn qiav it ('flow') and m qhish u ('<7high') distributions. The highest 
black curve shows the sum of the two, mq hoth ii ('fl>oth')- Notation is defined in the text. Blue: 
theory distribution ('TH'). Red: inconsistency cuts ('incons. cut') have been applied on m q i(\ ovl ) 
and m g ;(hi g h) (the resulting distribution is scaled up). Right: An inconsistency cut is applied on 
Tfiqii . The jet chosen for Tn qiaw ii is the one used in JTi qiow /(io W ) e t c - The inconsistency cut is seen to 
work very well for the m q i(\ ow ) distribution. 



The selection of the jet to go with the two leptons is more difficult. In the precuts at 
least three jets are required. If the squarks are considerably heavier than x!>, and the gluinos 
are not very much heavier than the squarks, then the two py-hardest jets are expected to 
come from the decay of the two squarks present in nearly all events. (More specifically we 
need m? — m_ n somewhat larger than m~ — m?) Along the entire SPS la line this is the 
case. The third jet is then expected to come from a gluino decay. Specifically, both for (a) 
and {(5) the correct quark is one of the two px hardest in 94% of the cases at the parton 
level. It is therefore reasonable to consider only these two jets. 

In a realistic setting all possible combinations will be investigated, along with all kinds 
of precuts. However, whatever the mass situation, one will not be able to select the correct 
jet by some simple cut. Combinatorial background from the jet selection procedure is 
something one will have to live with. 

In Fig. ^ (left) the different-flavour-subtracted m q u distribution at (a) is shown for 
various jet selections. The two lower black curves are obtained by selecting from the two 
PT-hardest jets the one which gives the smaller (solid) and larger (dashed) value of m q u, 
for simplicity referred to as m qiaw u ('flow') and rn qhighU ('ghigh')i respectively. The highest 
black curve shows the sum of the two, m qboth u ('fl, th')- In blue and marked 'TH', with 
the same normalisation as Tn qiaw u, is the parton level theory distribution, as generated by 
PYTHIA so that it contains widths, production biases etc. 

It is possible to get a reasonable first estimate of the endpoint both from rn qiaw ii and 
Tn qhoth u. The shape of the theory distribution is however not reproduced by any of the 
three experimental distributions (in black), which is quite natural since they all contain 
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both combinatorial and lepton-correlated background, and two of them are also biased by 
the jet selection procedure. For (a) the lepton-correlated background comes mainly from 
sleptonically decaying Xrj's that do not have a qL or b as a parent. As can be seen from 
Table |3] together with Fig. ||, \\ comes from qi or b in 78% of the times. If, for a signal event, 
we assume 50% chance of picking the correct jet, a different-flavour-subtracted distribution 
will consist of 39% signal, 39% combinatorics and 22% lepton-correlated background. For 
(/3) the fractions are 35%, 35% and 30%. For these numbers X2 1S assumed to be the only 
source of lepton-correlated background. 

We have investigated methods to remove or model this combinatorial background. For 
example, the sample may be purified [25] by using events where only one of the jets can 
reasonably be assigned to the signal quark. Consider the case where a measurement of 
m^ ax has been made, and one wants to now measure m^a^y In constructing the m^( low ) 
distribution we must choose one of the two highest px jets, and with no other input will 
choose the wrong jet 50% of the time, giving a combinatorial background. However, in 
some events, one of the jets will be such that its intepretation as the signal quark would 
contradict the previous measurement of rn^ x , i.e. using it to construct m q u would give 
m qii > m ^if x - m these events, we can be reassured that in choosing the other jet we 
have chosen correctly. If neither of the jets in the event gives m q u in contradiction with 
the previous measurement, then we cannot be certain that we have the correct jet, and 
therefore discard the event. In other words, we only use events where only one choice of 
signal quark is consistent with other endpoint measurements. This allows one to build a 
purified sample where for signal events the wrong jet is only chosen if the correct jet is 
not one of the two £>T-hardest jets. There will however still be background events in the 
sample. 

The red curve of Fig. |l^ (left) is constructed from the m q u distribution by means 
of 'inconsistency cuts' on m q ia ovr ) and m-^high) where we insist that only one of the two 
highest px jets has m^j^ < 310 GeV and Wg^high) < 390 GeV. Approximately 1/4 of 
the events pass these cuts; the resulting curve is normalised to have the same number of 
events as the m gio u distribution. The similarity with the blue theory curve has improved, 
but is still not excellent. In particular there is a large tail towards higher values, which in 
practice makes it less useful for accurate endpoint determination. 



In this respect a better result is obtained for the m g /(i ow ) distribution, shown in Fig. |12 



(right) with an inconsistancy cut requiring only one jet with m q u < 440 GeV, and normal- 
isation as rn qiow i(i ow y The red curve matches quite well with the theory curve, and there 
is no disturbing tail. In fact the red curve is the one we later use to estimate the endpoint 
for this distribution. 

While it is certainly reassuring to regain the theoretical distribution, it is still not clear 
what is the best way to obtain an accurate measurement of the endpoints. One would prob- 
ably combine many different methods, looking for convergence as well as inconsistencies. 
For instance, even though we should choose not to fit the distribution from the inconsis- 
tency cut, since it has less statistics, it teaches us that a linear fit to m qiow i^ ow ^ from 150 
to 300 GeV, which might otherwise seem reasonable, would not at all be a fit to the signal 
since it is close to linear only in the region 240-300 GeV. 
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Another approach to obtaining a distribution closer to the original one relies not on 
purification of the actual sample, but on finding an estimate for the background, combi- 
natorial and/or lepton-cor related, and then subtracting it from the original distribution. 
Statistics may be better preserved in this way, but systematics are introduced. Some at- 
tempts have been made with 'mixed events', i.e. the combination of the lepton pair from 
one event with jets from other events. The idea is that since the lepton and the jet sectors 
of mixed events are necessarily uncorrelated, they may mimic both types of background, 
lepton-correlated and combinatorial, where also the lepton pair is only weakly correlated 
with the jet. 

A first complication of this method is encountered when the four-vectors of different 
events are to be combined into mixed events. What rest frame is the appropriate one? 
This already points to the inexact nature of such a method. Its performance is however 
promising, although somewhat variable, as will be demonstrated later in this section. 

The mixed event sample is constructed from the events which make up the distri- 



butions, i.e. those which pass the precuts listed in Sect. p7T . An additional requirement 
of having exactly two opposite-sign leptons is imposed to avoid ambiguities in the lepton 
sector. No major differences are however expected from leaving out this constraint. Each 
event then provides two leptons, uniquely defined because of the additional cut, and its 
two pT-hardest jets. The four-vectors of the selected leptons and jets are here taken in the 
rest frame of the lepton pair. Due to the prominent role of the lepton pair, this choice 
is not unreasonable. On the other hand it is not unique. Both the laboratory frame, 
the rest frame of the entire event or the rest frame of the leptons together with the two 
pr-hardest jets are viable candidates. A dedicated study of the characteristics of these 
different possibilities would be worthwhile. 

The mixed sample is then constructed by consecutively combining the lepton sector 
of one event with the jet sector of another. In our study each lepton sector was combined 
with the jet sector of five other events. By increasing this number, very high statistics can 
be obtained for the mixed sample. For each combination of events all relevant masses were 
constructed, m giow /z, m qh . gh u, m qbothll , m giow ,(i ow) etc. 

Mixed samples for same-flavour and different-flavour events are constructed separately. 
In the end their difference is taken, in line with the previously described different-flavour 
subtraction. The resulting mixed distributions are the ones used later in this section. 

5.4 Multiple squark masses 



The theory distributions shown in Fig. 10 are only for one squark mass per scenario. In 



reality four different squark masses, mj^/m^, m,u L /mc L , and will contribute 
to the distribution. The blue theory distribution of Fig. 12 is therefore the sum of four 
separate distributions, each with a different endpoint, different normalisation and similar 
though not identical shapes. 

While jets from the three lightest quarks, and for most purposes also from the c-quark, 
will be indistinguishable in ATLAS, jets from 6-quarks will be identified with a certain 
probability. The expected rejection factors for incorrectly identifying jets from lighter 
quarks and gluons as 6-jets, are for high luminosity operation given in the ATLAS TDR 
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|2C], Fig. 10-41. Low/high 6-tagging efficiencies come with high/low rejection factors and 
allow for high-purity ft/non-fr-samples, respectively. The higher the purity, the smaller the 
sample. In this analysis we have used the following simplistic 6-tagging prescription for 
both purposes: For a 6-tagging efficiency of 50% the rejection factors against jets from 
gluons/three lightest jets and from c-jets are set to 100 and 10, respectively. 

With 6-tagging one can to a certain extent separate the qi and the b distributions, 
thus opening for a disentanglement of the squark masses. Nevertheless, even though a high 
purity separation has been accomplished, each of the two distributions will still contain 
contributions from two squark masses. Typically a kink can be observed at the position of 
the lowest endpoint. For mSUGRA scenarios di/si and ul/cl only differ by a few GeV, 
so the kink will appear very near the end of the distribution. Since the proton contains 
more u than d-quarks, ul will be produced at a higher rate than the heavier di- This 
reduces further the visibility of the kink. Then with the general smearing due to physics 
and detector effects in addition to background near the endpoint, it may be very difficult 
to identify such a kink. In case of the two 6-squarks the separation will be larger. Whether 
it is possible to identify it or not depends on the rate of b production as well as the level 
of impurity from g^-events in the 6-tagged distribution. 

5.5 Invariant mass distributions 

SUSY processes for (a) and (/?) as well as the Standard Model background have been 
produced for 300 fb" 1 . This corresponds to 3 years at design (high) luminosity. The mass 
distributions of the available edges for SPS la (a) and (/?) are shown in Figs. |l3|-14. 



For all plots the black points with errors show the total different-flavour-subtracted 
distribution ('SF-DF'). Solid green marks the SUSY background ('SUSY'), and in green 
with error bars the Standard Model background ('SM') is shown. The solid blue curve then 
shows the original theory distribution ('TH'), normalised to the different-flavour-subtracted 
distribution. The fitted function appears in red. In cases where mixed events are used to 
model the background, the smooth function fitted to the high-statistics mixed-event sample 
is shown in dashed red. When additional distributions are plotted, they will be described 
in the accompanying discussion. 

For each distribution the endpoint estimation will be discussed. In most cases the 
edges are fitted to a straight line in combination with a simple background hypothesis, 
and in some cases convoluted with a Gaussian distribution. At (a) this procedure gives 
numbers in reasonably good agreement with the nominal values. At (/?), where the SUSY 
cross-section is much smaller and also the branching fraction of the signal is reduced, 
the estimated endpoint values depend more strongly on the fitting method chosen. To 
control and reduce this systematic effect, a better understanding of the whole chain is 
required; physics effects, detector effects, multiple masses at different rates, background, 
precuts. After some years of LHC operation one can expect these issues to be understood 
sufficiently that the systematics of endpoint estimation is controlled and corrected for in 
the fitting procedure, up to some small uncertainty. 

If this is achieved, it is the statistical error of the endpoint values, in combination 
with the uncertainty on the absolute energy scale, expected to be 1% for jets and 0.1% for 
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leptons, which determine the precision at which masses can be obtained for the different 
scenarios. 

In the discussion of the distributions we therefore have the main focus on the statis- 
tical uncertainty of the endpoint determination, but are also concerned with the present 
magnitude of the systematic uncertainty. For some of the distributions it may seem bold 
to state that the systematic fit uncertainty will be reduced far below the statistical one. 
One shall however remember that once such distributions become available, a great effort 
will go into investigating them. 

SPS la (a) 



mu\ The mu distribution at (a) is shown in Fig. 13 (top left). It is very close to the 
expected triangular shape. To the right of the endpoint the Z peak is clearly visible. These 
Z's come from SUSY decays. The bump at low mu is discussed in Sect. |5T 



The distribution is fitted with a straight line convoluted with a Gaussian, 



\Am u , m u <mf ax 1 f°° . hiXi) 2 

9{niu) = < maY , /(mH) = -=-/ dm u g(m u )e ^ (5.1) 

[0, mu > m™ ax V27TO- Jo 

Both the endpoint, the width of the Gaussian and the normalisation are kept free in the fit, 
resulting in an endpoint value of 76.72 GeV at a statistical error of 0.04 GeV. Compared 
to the nominal value this is 0.35 GeV too low. When a more appropriate model for loss of 
lepton energy is implemented in the fit function, such a systematic shift is expected to be 
reduced to the level of the statistical error or below. The fit results are quite stable with 
respect to variations in bin size and range of fit. We thus arrive at the well-known result 
that the mu endpoint is expected to be estimated to high precision. 

m q ii: As for all cases where jets are involved, the m q u distribution comes in many ver- 
sions, depending on the jet selection procedure as well as the use of 6-tagging informa- 
tion. Plotted in Fig. [l3| (top right) is m qiow ii (see Sect. |5.3| ) with no 6- veto. Mixed events 
have been used in an attempt to model the background. A 6 th degree polynomial was 
first fitted to the high statistics mixed event sample. Then, with the normalisation kept 
free, the polynomial was combined with a straight line for the signal part. This proce- 
dure returns m^f* = 427.7(0.9) GeV, statistical error in parenthesis. The nominal value is 
425.9/431.1 GeV for ul/cIl- Alterations to the bin size and fit range has a 1-2 GeV effect 
on the fit value, but the statistical error remains the same. 

The dashed red line shows the shape of the mixed event sample. It matches the m qiow u 
distribution very well at values beyond the endpoint, justifying its use, and is sufficient to 
obtain a reasonable endpoint fit. Similar results can be found by modeling the background 
tail with an exponential or a polynomial. If, at some later stage, the entire distribution is 
to be compared with the theory distribution, the size of the background must be known 
for any m„u value. An exponential or a polynomial based on the upper tail can never be a 
good background model for values somewhat below the edge. It is here the mixed events 
may have a role to play. In the particular case of the m q , u distribution, the mixed sample 
actually appears to be a good background estimate at all values. For m qhath u and some 
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Figure 13: Invariant mass distributions for SPS la (a). See the text for details. 



other distributions, it does not do so well. Nevertheless, there are many ways to construct 
the mixed sample. Further study of the method seems worthwhile. 
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It is tempting to suggest that it is mainly the ul/cl endpoint we measure, and that 
the extra events just above the endpoint (see the figure), make up the hard-to-detect 
kink from the cIl/sl edge. At the present level of detail it is impossible to say if such 
an effect could be isolated to give an additional measurement for the heaviest mass. It 
does seem difficult, though, as the kink will be washed out by other effects, e.g. sparticle 
widths. At (a) the intrinsic width of ul is 5.3 GeV, see Fig. ^| (right), to be compared 
with mj^ — rriu L = 5.8 GeV. Also detector effects will result in a general smoothing of the 
distributions. 

m qi(iow) : For m q in ow -\ the mixed sample is fine above the endpoint, but overestimates 
for lower masses. While the edge can be fitted and the endpoint measured by a mixed 
sample subtraction, due to its good behaviour in the edge region, we have here instead 
used the inconsistency cut m qh . h u > 440 GeV to purify the sample, Fig. [l3| (middle left). 
The ability of the inconsistency cut to bring the distribution very close to the original one, 
was already discussed in relation to Fig. [T^ . 

For a zero background hypothesis a straight line fit gives w^j^ = 300.7(0.9) GeV, to 
compare with the nominal 298.5/302.1 GeV for ui/di- If the few bins around the endpoint 
are also included, the value increases by 1-2 GeV. Whether these high-mass events are 
signal or background is not easy to tell from the given distribution, since there is virtually 
no background structure to extrapolate from. This is because a 'consistency cut' has 
been imposed which requires rn q . u < 440 GeV, in accordance with the already measured 
m q ii endpoint. From Eq. ( |4 . 1 0| ) this implies rn q i^ ow ^ < 440/y / 2 GeV = 311.1 GeV. The 
consistency cut takes away a large part of the background, but also has the effect that it 
becomes difficult to see what structure the m q i(i ow } background has. If the consistency cut 
is dropped, the usual background tail appears and can be modeled, although at the cost 
of a slight increase in the statistical error. 

m qi(high) : Following the same procedure as for m q u the background of the m g z(high) dis- 
tribution was modeled by the mixed event sample. In Fig. [l3| (middle right) the relevant 
distributions are shown. Also the result from subtracting the mixed sample is shown 
('mixed subtr.'); the lower black points with error bars. This subtracted distribution fol- 
lows the original theory distribution (blue) closely in the edge region, but overestimates at 
lower values. 

In the range m g z( high ) € (320, 550) GeV the endpoint was estimated to 374.0(2.0) GeV. 
The nominal value is 375.8/380.3 GeV for ul/<1l. At (a) the m^high) edge consists of 
two parts. This is clear from the theory distribution (blue), but also in the reconstructed 
distributions (black). In the fit only the lowest near- linear stretch was used. It is clearly 
incorrect to apply a straight line for the whole edge, but if done, the statistical error would 
be reduced to ~ 1 GeV. When a good signal function is at hand, the whole edge will be 
described. It may therefore be reasonable to expect a statistical error of 1 GeV rather than 
2 GeV. 

m qii(e>-) : The r ra q n^Q > iL^ distribution differs from the previous ones in that a minimum 



is to be measured. This has two important consequences. One can be seen from Fig. 11 
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(bottom left), even though this plot is for the unconstrained m q u distribution. The lepton- 
uncorrelated background, modeled by the different-flavour distribution (dashed red), sits 
mostly at low values, far away from any maximum edge, but quite near the m q n^ e> ^^ 
threshold. While this background type is removed (statistically) by the different-flavour 
subtraction, the price is an increased statistical uncertainty which directly will affect the 
statistical precision of the fn^^Q > ^ estimation. 

The other consequence of measuring a minimum is related to the multiple squark 
masses. As noted earlier, the distributions can be purified using 6-tagging information. For 
the measurements discussed so far we have considered the total sample. This is because 
the gx-squarks are so much heavier than the 6-squarks that the upper edge regions are 
practically free of 6-events anyway. A b- veto has therefore no effect, at least not at our 
level of precision. For the threshold measurement it is the other way around. Here the 
threshold related to the lightest squark, b\, will be the threshold for the total distribution. 
A 6-veto may therefore be important in order to obtain the threshold for the non-6-events. 

A third complication is that the edge is very non-linear, as can be seen from the theory 
distributions in Fig. 10. Also, the fact that this non-linearity is of the concave type, could 
make it very difficult to see where the edge ends and where smearing from various sources 
takes over. This is all the more relevant since the 6-events come at a significantly smaller 
rate than g^-events (ratio ~ 25%), and may therefore appear as just another contribution 
to the lower tail, not easily distinguishable from other effects. 

The J7i 9 zz(0>i) distribution is shown in Fig. |l3| (bottom left). The upper distribution 
('SF-DF') shows the non-6-tagged m 9 zz(0>2t) sample. The lower black points ('SF-DF(fc)') 
mark the 6-tagged sample. Only approximately 50% of the events which contain a b 
are actually 6-tagged. This means that the non-6-tagged distribution contains a 'hidden' 
6-sample, statistically similar to the 'visible' 6-sample both in size and shape. (There 
would be a shift to the left, though, since in the reconstruction 6-tagged jets have larger 
recalibration factors. This is not accounted for here.) Analogous to the different-flavour 
subtraction it is therefore possible to get the non-fo-tagged distribution closer to the original 
(/-distribution by subtracting the visible 6-sample, possibly scaled up or down according 
to the 6-tagging efficiency applied. A slightly different approach would be to use a much 
harder 6-veto and no subtraction. The price would be a significant reduction in sample 
size. 

However, the main problem in the threshold region is not the 6-contamination, but the 
lepton-correlated SUSY background shown in solid green. It comes primarily from events 
where x% decays sleptonically but does not descend from qi or b. It is this background 
which needs appropriate modeling. The difficulty with having a concave theory distribution 
together with an unknown background is evident from the figure. It is really difficult to 
separate the background from the signal. 

As a simple estimate of the statistical uncertainty of the edge position, a straight 
line fit somewhat away from the threshold was performed on the non-subtracted sample 
('SF-DF'), returning a statistical error on the endpoint around 1.8 GeV. If the subtracted 
sample (subtracting 'SF-DF(6)') is used, the error increases. The actual endpoint obtained 
from such a simple zero-background hypothesis is of course incorrect. At present it seems a 
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bit optimistic to expect the systematics of the fit to be dominated by the statistical error. 

m bii(0>^) : The w<Mi(6i>|) distributions are shown in Fig. (bottom right). The different- 
flavour-subtracted distribution ('SF-DF') is shown in black with error bars. The same- 
flavour ('SF') and different-flavour ('DF') curves are also plotted. In the threshold re- 
gion the same-flavour and different-flavour distributions are considerably larger than their 
difference, which necessarily gives large statistical uncertainties in the different-flavour- 
subtracted samples. 

As for m q n(Q > nj more studies are required to control and reduce the systematic error 
induced by the background and the nonlinear theory curve. A linear fit was performed in 
the near-linear region to measure the precision with which the position of the edge can be 
found. Variations in the bin size and the fit range yield systematic errors within 1 GeV 
and statistical error between 4 and 5 GeV. 



rribu, m bi (i ow ), m W ( high ): Given that the background can be controlled and a good 
signal function can be found, rribu^e^) is a good distribution since the lowest threshold is 
given by &i, which also has a considerable rate. For the edges mja, mbi(\ow) arm f&w(high) one 
problem is that it is 62 which defines the outer endpoint for the 6-distributions. Since this 
squark only contributes ~ 22% of the 6-events, the edge above the b\ endpoint may look 
like background. (However, since we expect the heavier 6-squark to come at a considerably 
smaller rate as it has less phase space for production, and also a smaller branching ratio to 
the wino-like x^j due to its larger right-handed component, we can explicitly look for such 
a small tail-like edge.) 

A more serious problem is the background from other SUSY events. Typically these 
have the decay Eq. ( [3.3j ) in one chain, providing two leptons, and g — ► bb or g — ► t\t in 
the other, providing 6-jets. The mass distributions of this background stretch sufficiently 
beyond the 61 endpoints that the resulting edge structures become difficult to analyse. The 
positions of these three 6-edges can typically be measured with a statistical precision of 
3-4 GeV, but it remains unclear whether the systematics of the edge can be sorted out. 
Another approach could be to require a lower 6-tagging efficiency for these distributions in 
order to get higher rejection factors and purer, although smaller, samples. 



SPS la {ft) 

The mass distributions of (/?) are shown in Fig. |14|. Throughout, the total samples have 
been used. No 6-veto has been applied. 



ma: The mu distribution is shown in Fig. 14 (top left). It has the triangular shape with 
a Z-peak on top. Some datapoints in the Z-peak lie above the range of this plot, and can 
be viewed in Fig. [IT] (top right). 

As was also the case for (a) there is a bump at lower values. The distribution above the 
Z-peak was fitted with the same Gaussian-convoluted triangular shape as earlier, Eq. (|5.1|), 
giving m™ ax = 137.4(0.5) GeV which is 0.5 GeV lower than the nominal value. 
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Figure 14: Invariant mass distributions for SPS la ((3). See the text for details. 



m q u: The Tn qiow ii distribution is shown in Fig. [l4| (top right). With the mixed sample 
as a rough background estimate a straight-line fit gives an endpoint at 640 GeV with a 
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systematic shift of 1 GeV from varying the bin size and fit range. The statistical error is 
5 GeV. If, however, the same fitting procedure is applied to the m qboth u distribution, the 
endpoint value increases to 655 ± 2 GeV and the statistical error is 9 GeV. 

Surely more study would bring these values closer, and optimally have them converge 
near the nominal value of 649.1/652.5 GeV for ui/di- For later use we take an optimistic 
statistical error of 5 GeV, but also include a systematic fit error of 3 GeV for a more 
conservative estimate. 

m qi(iow) : The m,j low i(iow) distribution is shown in Fig. [L4| (middle left). Both a consistency 
cut m qiov/ u < 670 GeV and an inconsistency cut m qhigh u > 670 GeV are used. A straight- 
line fit with no background hypothesis gives a statistical error around 6.3 GeV. The actual 
fit value is 443 GeV, which overshoots by a few GeV since the background plateau has not 
been included. The nominal value is 436.6/438.9 GeV for u^jd^. 

m qt(high) : The m giow ^hi g h) distribution with the consistency and the inconsistency cut is 
shown in Fig. [li] (middle right). The mixed sample was again used to roughly model the 
background under the edge. A straight-line fit gives an average value of 520.5 GeV with 
systematics from binning and fit range of 3 GeV. The statistical error is at 5.5 GeV. 

The nominal value is 529.9/532.7 GeV for u^/dL, some 10 GeV above our estimate. 
One could argue that the current endpoint measurement has considerable uncertainties 
and that this discrepancy is not dramatic at the present level of detail. However, such an 



underestimation is actually to be expected. In Fig. 1C the theoretical ^/(high) distribution 
is shown for ul at (/?). There is a long vertical fall towards 517 GeV (for di, it is at 
519 GeV), then just before the bottom is reached, a small foot appears, as anticipated in 
Sect. |4.1| , and takes us up by 11-13 GeV. To detect such a small foot would require more 
statistics than is available at (/?). Experimentally it is therefore expected to get an endpoint 
near 517-519 GeV. The incorrect endpoint measurement will have important consequences 
for the determination of masses from the endpoints. This situation is further discussed in 
Section ||. 

m qii(0>^) : The m q ii^ d> ^ distribution is shown in Fig. |l4| (bottom). The same-flavour 
('SF') and different-flavour ('DF') distributions are shown in solid and dashed black. 
Clearly, there are large uncertainties from the different-flavour subtraction. The statistical 
uncertainty of the different-flavour-subtracted sample was estimated with a straight-line 
fit, giving 13 GeV. In addition there will be a systematic error, here conservatively set to 
10 GeV. 

Endpoint measurement values 

The results of the endpoint estimation for (a) and (f3) are summarized in Table ||. The 
last column contains an estimate of the systematic error from different fitting techniques, 
ranges and bin widths. These values are not used in the following, but are included for 
completeness. 

The column with the heading 'Energy Scale Error' shows the expected error on each 
endpoint estimation from the uncertainty on the absolute energy scale for jets and leptons. 
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Table 4: Endpoint values found from fitting the edges in Figs. [TH-fl, for 300 ft." 1 . The nominal 
values correspond to the mass of ul, which due to the proton content is produced at higher rates 
than the heavier <Il- For the thresholds no fit values are shown, only the errors. This reflects the 
fact that a reasonable fit function is lacking for this edge. 



This effect has not been taken into account in the simulation. The uncertainties of the 
energy scales are here set to 1% for jets and 0.1% for electrons and muons, see Ch. 12 of 



|2C]. For an invariant mass which consists of only jets or only leptons, this will give the 
same uncertainties, 1% and 0.1%, respectively. If the invariant mass is constructed from 
one jet and one lepton, the endpoint uncertainty is 



<r{m q i) _ (r(m 2 gl 



rriqi 2m? 



ql 




0.50% (5.2) 



where Ej and E\ are the jet and the lepton energies, respectively. For an invariant mass 
involving a higher number of jets and leptons, the error on the endpoint value from the 
energy scale uncertainty will be different for each event. The error of m q u will depend on 
the relative size of the three terms on the right-hand side of Eq. ( [4.10 ). Since at (a) and ((3) 
we are in the region where the mass ratio rriq L /m^o dominates the two other mass ratios, 
see Eq. the quark will usually be very energetic, leaving one or both m q i terms 

to dominate. This is particularly true at large values, so near the edge of m q u one can 
show that the energy scale error will result in an endpoint error between 0.35% and 0.5% 
for each event. For ra q n^Q > jL) the average energy scale error will be slightly lower in our 
two scenarios. For simplicity we have set the energy scale error to 0.5% for all endpoints 
involving jets. 
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6. Extraction of masses from edges 



6.1 10,000 ATLAS experiments 

In the simulation study described in Sect. || values for the endpoints and their statistical 
uncertainties were found, together with a 'systematic fit uncertainty'. Although not so far 
from the nominal endpoint values, the fit values in Table |I] are somewhat uncertain due to 
the as yet not-understood systematics of the fitting procedures. Also, the systematic error 
on the energy scale has not been addressed. 

Assuming that one will eventually be able to control the systematics of the fitting, only 
the statistical errors together with the systematic error from the energy scale uncertainty 
will be what determine the LHC potential to measure the SUSY masses. To estimate this 
potential, consider an ensemble of typical LHC experiments, i.e. where the deviation of 
each endpoint measurement from the nominal value is based on a Gaussian distribution 
of width equal to the statistical error estimated for that endpoint, as well as a jet/lepton 
energy scale error picked from a Gaussian distribution for each experiment, in line with 
what is done in [31], 

^exp = gnom + ^stat + ^scale (g 

Here Ei denotes the position of the z th endpoint. The coefficients A and B are picked from 
a Gaussian distribution of mean and variance 1. Each experiment will pick as many A's 
as there are endpoint measurements as well as one B for the ma endpoint and one for the 
endpoints involving jets, thus neglecting the effect of the lepton energy scale error on the 
latter. 

When a set of edges E exp has been found, the task is to find the masses m which 
best correspond to the measurements. If only four endpoints are measured, the inversion 
formulae straight away return the possible mass combinations. If more endpoints are 
available, no mass combination will in general reproduce the edge measurements, and a 
numerical approach is required, where the measurements are weighted according to their 
uncertainties. Note that in this procedure we do not make use of the fit values given in 
Table | 

It should also be emphasised that the systematics of the endpoint measurements are 
here assumed to be under control, i.e. the 'Syst. Fit Error' of Table |] is neglected. The 
precision we will find in this section and the next for the determination of masses and mass 
differences at the LHC must be understood in this context. If the endpoint systematics 
turn out to be comparable to the combined statistical and energy scale errors, then the 
precision will be worse. 

6.2 Mass estimation via £ 

In our case, where the jet energy scale error produces a correlation between the endpoint 
measurements, the method of least squares is appropriate. The best mass estimate rh is 
then the one which minimises the function 

£ = [E ex P - E th (m)] T W[E ex P - E th (m)] (6.2) 
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where E th (m) contains the theoretical edge values for a set of masses m. The weight matrix 
W is the inverse of the correlation matrix or error matrix of the observations, which is given 
by the variances and covariances of the endpoint measurements, 

(W- 1 )^ = of| at + a^ c = (of at ) 2 + (crf al0 ) 2 

(W~% = cr|f lc = (£f p £f p ) - (Ef p )(Ef p ) = of al Vf ale , % + j (6.3) 
(W- 1 )*! = 0, » ^ 1 

where j = 1 refers to m™ ax , which to a good approximation is uncorrelated with the other 
measurements. The covariances are similar to the variances in size, and so cannot be 
neglected. If the endpoint measurements were uncorrelated, W would become diagonal, 
and the least-squares method would reduce to the normal x 2 minimum method. 

The ensemble distributions obtained by such a procedure can be interpreted as prob- 
ability density functions. Prom these the 'inverse probability problem' can be addressed, 
which is that of stating something about the true masses on the basis of the ones obtained 
in one experiment. We will be interested in the mean values of the ensemble distribution, 
their standard deviations, skewness, as well as the correlation between masses. 

6.3 Minima of £ 

Because many endpoints are given by different expressions for different mass regions, see 
Eqs (|4.4j)-( |4~5| ), the minimisation function £ is a composite function. For the endpoint 
measurements used in this paper, £ is made up of nine individual functions, £(ij), one 
for each of the nine regions Considered separately each function E^j) has one or 

more minima. For these to also be minima of the composite function ('physical minima'), 
they need to be situated in the region of validity ('home region') of the corresponding 
function. Physical minima can also occur on the borders between regions, in which case 
they will be referred to as 'border minima'. 

If the threshold endpoints are left out, there are four measurements for four masses. 
The clear failure of the endpoint measurements of SPS la (a) and ((5) to comply with 



Eq. ( 4.11 ) already discards the three regions where these four measurements are not suffi- 
cient to determine the masses. In each of the other six regions the minima can be sought 
by use of the inversion formulae. Such solutions correspond to E = 0. In cases where no 
physical solutions are found in this way, border minima exist at £ > 0, and will have to be 
found by a least square minimisation. 

When the threshold measurement is included, the system of equations becomes over- 
constrained. This will give a non-uniform increase in the value of £, which may destroy 
or create minima. Another effect will be to move the minima of £ around in mass space, 
possibly moving them into or out of their home regions. One way to picture the effect is 
to 'tune in' the new measurement by letting its uncertainty go from infinity, in which case 
the measurement has no effect, to the value specified in Table ||. The masses and height 
of each E minimum will then move continuously from the old to the new position. 

Even though composite, E is continuous, so its realisations in two neighbouring regions 
attain the same value at their common border. Assume that the endpoint measurements 
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are such that one of the realisations has a minimum at the border (not a so-called border 
minimum). Consider first the case of no threshold measurement. Since E = 0, also the other 
realisation must have a minimum at the border. If now one of the endpoint measurements 
is shifted up or down, the two minima will be driven off the border and also separated in 
mass space. To which side of the border they move, and whether they go to the same side 
or not, will depend on the actual parameters. If the minima are on the same side, as is the 
case at (/?), only one of them will be a physical minimum. If they are on opposite sides, as 
is the case at (a), either both will be physical minima, or neither, in which case there will 
be a border minimum. 

This means that if a mass scenario is situated close to one of the borders, the endpoints 
it produces may also have been produced by a set of masses from the neighbouring region, 
provided the minimum there is a physical minimum. 

If the threshold measurement is added, this picture is no longer exact. Since £ does not 
vanish at the minimum, the two regions will in general no longer have a common minimum 
at the border. However, the threshold measurement often has less weight than the other 
measurements, so the above picture still has some validity: Near a border two minima will 
be lurking. Both, one or none may be physical minima. 

6.4 SPS la (a) 

If we neglect the threshold measurements to start with, and thereby also wait with the b\ 
mass, there are for each 'experiment' two solutions. One solution is in region (1,1), which 
is the home region of the nominal masses, and one is in region (1,2). This is an example 
of the situation described above: SPS la (a) has 2m? /(m 2 ~ Q + m? ) = 1.01. If this ratio 
becomes less than one, the region changes, as is seen from Eq. ( |4.5| ). A reduction of mj R 
by 0.7 GeV would put the mass set on the border to region (1,2). In both regions the 
mass distributions are close to Gaussian. The ensemble means of the home region solution 
essentially equal the nominal values. The (1,2) solution has central values some 15-20 GeV 
below the nominal ones. Without additional information both solutions are equally good 
and it would not be possible to determine which one to choose. One would have to state 
that the SUSY masses are summarized by either the (1,1) set or by the (1,2) set. 

If the non-6 threshold endpoint is included as a fifth measurement in a least square 
minimisation, the situation changes. While a (1,1) minimum exists for practically all the 
'experiments', the occurence of a (1,2) minimum is now slightly reduced to 85-90%, but 
it is in the £ value of the minima the effect is most apparent. The overconstraining fifth 
measurement lifts the two minima asymmetrically from zero. As could be expected, the 
(1,1) minimum is more often in accordance with the data, but there is always a non- 
negligible (1,2) contribution. The difference in E value between a given minimum and the 
global minimum, AS, is a measure of the relevance of the minimum. Table || shows how 
many minima are available, on average, for various AS cuts, and how these are shared 
between regions (1,1) and (1,2). 

If only the lowest minimum is chosen, the wrong solution is returned in 10% of the 
experiments. However, if two minima exist and are close in S value, one would have 
to consider both. In a certain fraction of experiments, depending on the AS cut, there 
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would thus be two solution sets, e.g. for AS < 1 we will have two solutions in 12% of 
the experiments. Whether or not it would be possible to select one of the solutions, and 
preferably the correct one, hinges on other measurements. In this case, where the masses of 
the two sets are quite close, they might be very difficult to distinguish, by e.g. cross-section 
considerations. 
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Table 6: SPS la (a): Minima for AE < 1 in regions (1,1) and (1,2). Ensemble means, (to), 
and root-mean-square distances from the mean, cr, are in GeV. The three lightest masses are 
very correlated. The mass of qL is fairly correlated to the lighter masses, but mg is essentially 
uncorrelated. The distributions are very close to symmetric. 

The upper part of Table |6| shows the en- 
semble means of the masses, (m), the root- 
mean-square distances from the mean, a, and 
skewness 71 4 of the two solutions for AS < 1. 
The values are relatively stable with respect 
to the AS cut: The same table for AS < 99 
would for (1,1) show a decrease in the masses Table 5: Number of minima for various AS 
by 0.1-0.2 GeV, and for (1,2) an increase by cuts and their whereabouts. 
1-1.3 GeV. 

The inclusion of the threshold measurement has very little effect on the ensemble 
values of the (1,1) solution. For the (1,2) solution, to better comply with the additional 
measurement, the masses have increased, and are now 10-15 GeV below the nominal ones. 
Also the rribii threshold was included in the fit which returned the values of Table ^. It is 
measured with much less precision than the other endpoints, so its inclusion has practically 
no effect on the other masses, only on , for which it is the only measurement here. 

The fact that the ensemble means of (1,1) reproduce the nominal values, relates to the 
good average performance of the ensemble of experiments. The probability of doing well 
with only one experiment relies in addition on the spread of the ensemble values, given 
by a. For SPS la (a) the high precision of the endpoint measurements translates into rather 
small a values. From the table we see e.g. that in ~ 68% of the experiments the mass 
from the (1,1) solution will lie within 3.8 GeV of the nominal value. The root-mean-square 
4 Skewness is defined by 71 = ^/(a^) 3 ^ 2 = {13/a 3 , where im = (x — x) 1 is the i th moment about the 
mean. 
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Sparticle masses and mass differences [GeV] 

Figure 15: Sparticle masses and mass differences at SPS la (a) for solutions with AS < 1. The 
unfilled distributions in black show from left to right m^o, mj , m^o, mr and m*. for solutions in 
the nominal region (1,1). We will have such a solution in 94% of the experiments, see Tabic U The 
unfilled distributions in blue show the same masses for solutions in region (1,2). Such a solution 
occurs in 17% of the experiments, and the masses returned are lower. The smaller rate of the (1,2) 
solutions is reflected in the smaller area under the blue curves. The ratio of probabilities between 
(1,2) and (1,1) solutions is 17%/94% = 18%. The area under one of the blue curves is 18% of the 
area under the corresponding black curve. The filled distributions show from left to right mj R —m^o, 
m^o — m^o , — rn^o and wiq L — m^o . Again, the most populated distributions (black curves) 
are for solutions in region (1,1), the least populated (blue curves) for (1,2) solutions. For mass 
differences there is more overlap between the (1,1) and (1,2) solutions, in particular for mj^—m^o 
and m^o—m^o, of which only the lower parts of the distributions are visible. Mass differences are 
better determined than the masses themselves, reflected here by the narrower distributions of the 
former. The exception is which largely decouples from the other masses. 

distances from the ensemble means are in principle unknown, as seen from one experiment. 
They can however be approximated by the procedure of simulating 10 4 experiments, where 
the measured values play the role as 'nominal'. This will engender a systematic shift, but 
a and any skewness should be fairly well approximated. The root-mean-square distances 
from the mean values also have their counterparts in the la errors returned by the fit of 
each 'experiment'. To within a few percent they are found to be identical. This means that 
this information is available for the experiment actually performed. One can then make 
the inverse statement: For a given experiment one can with ~ 68% confidence state that 
the nominal value of m^o lies within 3.8 GeV of the mass returned. 

Due to the way masses enter in the endpoint expressions, the fit returns masses which 
have a strong positive correlation. If one mass is low at the minimum of the S function, 
so the others tend to be and by a similar amount. In the lower part of Table ^ ensemble 
mean and root-mean-square values of mass differences are shown. It is clear that the three 
lightest sparticles are very correlated. Fix one and the others are given very accurately. The 
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squark masses are less correlated. Also the results in region (1,2) are closer to (1,1) and the 



nominal ones when considering mass differences. Fig. 15 shows the ensemble distributions 
corresponding to Table ^. 

Because of this strong correlation between the masses, not only the mean values and 
their la uncertainties, but the entire error matrix should be considered if one wants to use 
the result obtained by this method as input for other analyses. 

A less involved solution would be to use less dependent variables, e.g. m^o to set the 
scale, then differences for the remaining masses, mr — m^o, m^o — m^o, m,fj T — m^o and 

° l R X\ ' X2 Xi ' qL Xi 

Due to the high cross-section, most of the endpoints are determined with high precision, 
which in turn gives narrow and approximately symmetric ensemble distributions. The 
masses are thus determined with quite high precision. As a result of the strong correlations 
between in particular the lighter masses, even better estimates can be obtained for other 
combinations of the variables, e.g. mass differences. At SPS la (a) there is however a 
fair chance that two sets of masses do equally well in the minimisation procedure. Other 
considerations must in that case be made in order to choose between them, or both must 
be kept. 

6.5 SPS la (f3) 

In combination with the theory plots of Section |^ we found in Section || that the m^m^^, 
value of {(3) would most probably on the average be underestimated by 11-13 GeV. We will 
find that this has a dramatic effect on the masses returned. However, for easier comparison 
with nominal values we first consider the situation without such a systematic effect. Then, 
afterwards, the impact of the mismeasurement will be shown, together with a way to mend 
the situation. 

As in the previous case we start out without the threshold measurement. Also here 
two solutions are available for all 'experiments', one in (1,1), the other in either (1,2) or 
(1,3). The nominal region for is (1,2), but it is quite close to (1,3). This becomes clear 
from an inspection of the 'border parameter', 



2 



Xi A 2 



(6.4) 



which is 1 on the border between these two regions. For {(3) we have b = 1.02. If mj is 
reduced by 2.5 GeV, this ratio becomes unity and the the mass set sits on the border, see 
Eq. (|4.5| ). The other point, (a), was near the border between (1,1) and (1,2). There, both 
solutions (or none) were available. Here, the derivatives of E are such that only one of the 
two solutions is available. In 71% of the cases we get (1,2), in 29% we get (1,3). While 
the 'low-mass' solution, (1,2) or (1,3), is in the vicinity of the nominal masses, the (1,1) 
solution, which is always present, usually sits at much higher masses, (m^o) = 514 GeV. 
Because the two solutions are so separated one may hope that the incorrect one will be 
sufficiently disfavoured by other measurements, e.g. cross-sections, that it can be discarded. 

For (a) the solution in the nominal region on the average reproduced the nominal 
values to within 0.2-0.5 GeV. Here the (1,2) j (1,3) solution has a mean of m^o at 183 GeV, 
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some 22 GeV above the nominal value. The most probable value of the ensemble distribu- 
tion is much closer to the nominal value. The distributions are infested with considerable 
skewness. On the way from endpoint measurement to mass determination a systematic 
effect which favours higher masses has been introduced. In statistical language our esti- 
mators of the true masses are not consistent: they do not converge to the nominal values. 
This of course has implications for the interpretation of the masses we obtain. What can 
be said about the true masses on the basis of the measured ones? We will return to the 
reasons for the skewness later. 

When the threshold measurement is included, the (1,1) minimum usually yields a large 
S value. Only in a small fraction of the experiments does it challenge the other minima. In 
the other sector there is either one minimum, positioned in (1,2), (1,3) or on the border (B), 
or there are two minima, in (1,2) and (1,3). These minima are usually in good agreement 
with the threshold measurement and have low S values. 
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5% 
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46% 


14% 


12% 


28% 


AS < 99 


2.3 


99% 


41% 


13% 


12% 


34% 



Table 7: SPS la (/3): Average number of minima and the fraction of experiments with the specified 
solution types, for different AS cuts. 

Table [7| shows the average number of minima for different S cuts. The three rightmost 
sections show the fraction of experiments which have the specified solution type. The two 
rightmost sections exclude one another. Either there is one solution in the low-mass sector, 
or there are two. For the one-solution case the whereabouts of the minimum is also shown. 
The home region of the nominal masses, (1,2), is seen to dominate. As it may well be 
possible to discard the (1,1) minimum on the basis of other observations, it is logically 
separated from the low-mass minima. E.g. in 13% of the cases, regardless of the low-mass 
solution type, there is a (1,1) minimum at AS < 3. To get the average number of minima 
shown in column 2, sum horizontally, adding twice the two-solution percentage. For small 
AS cuts the two rightmost sectors do not add up to 1. This simply means that in some 
cases the global minimum lies in (1,1), and no low-mass minimum is available in the given 
AS range. For the current set of endpoint measurements the (1,1) contamination is seen 
to be very moderate. However, the systematic fit error (column 6 of Table f|) is here 
assumed to be zero. If it should become impossible to obtain the threshold value with such 
optimistic precision, the fraction of (1,1) solutions at low S will grow rapidly. 

In Table || (m), a and 71 of the ensemble masses and mass differences are shown for 
the different solution types and the cut AS < 1. The masses of the (1,1) solution are much 
higher than what the low- mass minima give. Even though the distributions are broad, 
allowing for low values to occur, it is very rare that the masses stretch down to the low- 
mass sector. In section 4 of the table the low-mass one-solution values are shown. Since 
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Table 8: SPS la ((3): Nominal masses ('Nom') and AS < 1 ensemble distribution values for the 
three solution types. High-mass sector: The (1,1) solutions return masses far beyond the nominal 
values. Low-mass sector: For the one-solution case the values are based on the common distribution 
of (1,2), (1,3) and border (B) solutions. In the two-solution case the ensemble variables of both 
solutions are shown. Ensemble means, (m), and root-mean-square values, a, are in GeV. 



in such a case only one acceptable solution is available (discarding (1,1)), and since (/?) 
anyway is situated close to the border, it makes sense to show the combined distribution 
of the (1,2), (1,3) and border (B) minima. From Table ^ this situation is seen to occur in 
52%+18%+12%=82% of the experiments. The mean values of the masses lie some 15 GeV 
above the nominal ones. This is an improvement compared to the non-threshold situation, 
but it remains an undesirable feature. The mass distributions are skewed and the most 
probable value is found close to the nominal value. The root-mean-square values of the 
ensemble distributions are large, nearly an order of magnitude larger than at (a). 

The rightmost sections show the values for the two-solution type. The ensemble means 
of the two distributions do not differ too much, and they are much closer to the nominal 
values than is the case for the one-solution type. Since the values are rather close, it will 
probably be quite difficult to find other measurements which favours one of the sets. On 
the other hand, since the root-mean-square values are larger than the differences between 
the two solutions (also within one experiment), taking the average value, perhaps somehow 
weighted with the £ value, might be a possible compromise. 

Again, the mass correlation is very strong. This is evident from the lower part of 
Table ||, where the ensemble distributions of mass differences come with much smaller 
root-mean-square distance to the mean values than what the masses themselves do. They 
are also very close to the nominal values. Even more seems to be gained by using mass 
differences here than at (a). Fig. 16 shows the ensemble distributions for the masses of all 
solutions with AS < 1 and which lie in the regions (1,2), (1,3) or on their common border. 
See figure caption for details. 



Skewness 

The ensemble distributions are not symmetric. While the most probable values are close 
to the nominal values, the means lie above. For (a) the tendency of such an asymmetry is 
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Figure 16: Sparticle masses and mass differences at SPS la ((3). All masses of solutions with 
AE < 1 which lie in regions (1,2), (1,3) and on their common border are shown. From left to right 
the unfilled distributions show m^o , , m^o and rn,q L . The filled distributions show the narrower 
mass differences mr — m^o , m^o — m^a and m», — m^o . Skewness of mass distributions is visible. 



small, but for (/?) the effect is large. The reason why we naively would expect a symmetric 
distribution around the nominal masses in the first place, is that the endpoint measurements 
are generated symmetrically. For complex functions like Eqs. ( 4.13| )-( 4.3G| ) symmetric 
fluctuation of the endpoint arguments will produce near-symmetric variation of the function 
only for small fluctuations. As the arguments fluctuate more, the deviation from symmetry 
in the function values grows. At (a) the endpoint fluctuations are so small that the effect is 
negligible. For (/3), where the endpoint fluctuations are larger, the effect of the 'asymmetric 
propagation' is a noticeable increase of 3-4 GeV for the ensemble means. 

This is however not sufficient to explain the low-mass (ui^o) of 183 and 173 GeV 
(AS < 99) without and with the threshold measurement, respectively. 'Border effects' 
need to be considered. As described earlier, lies in (1,2) but close to the border to 
(1,3). First consider the situation without the threshold measurement. There is then 
always only one low-mass solution. If the (1,2) solution is physical, i.e. lies in (1,2), then 
the true minimum of ^n,s) also lies in region (1,2) and so is unphysical, and vice versa, 
as described in Sect. 6.3 . 

In Fig. 17 the mass of Xi 1S plotted as a function of the border parameter, 6, of 
Eq. (6.4), for both physical and unphysical minima of £(^,2) and Yln^\. The minima of 
are shown in red, from upper right to lower left. The ^(1,3) solutions are in blue. 
Filled boxes are physical solutions, i.e. ^(1,2) (red) for b > 1 and £(_/,#) (blue) for b < 1. 
Empty boxes are unphysical solutions. An asymmetry arises from the accidental fact that 
for both functions the lower masses tend to lie in the unphysical region. The average of 
the entire £(j,2) distribution, both physical and unphysical minima, returns 164 GeV, the 
nominal value plus the 3 GeV of the asymmetric propagation effect. It is then obvious that 
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Figure 17: Border effect between region (1,2) and (1,3); the mass of Xi as a function of the border 
parameter 6, see Eq. (6A), demonstrating border effects. Filled boxes represent physical solutions, 
empty boxes represent unphysical solutions. 



when the unphysical (1,2) solutions are replaced by the physical (1,3) solutions which lie 
at higher masses, the ensemble mean increases. Here, this effect brings the average value 
for m^o to 183 GeV, an additional increase of nearly 20 GeV. 

When the threshold measurement is added, the border effect is reduced, giving an 
average value of 173 GeV (for one solution). This is because a measurement will, unless 
there is any bias, on the average be conservative. It will try to keep the masses at their 
nominal values. Here, as the border effect is pushing the masses upwards, away from the 
nominal values, the threshold is holding back. 

If, in a realistic situation, a set of endpoints has been measured, and the resulting set 
of masses is found to lie close to a border, caution should be exercised. Ad hoc procedures 
seem necessary for such a case. At least one is in a position to be aware of the danger. It 
is probably appropriate to consider unphysical minima as well. 

Mismeasured w^mj gh ) 

As pointed out at the start of this subsection it is very likely that, for Tn^^g^ will 
be understimated by 11-13 GeV. Without the threshold measurement, the effect of this as 
returned by the inversion formulae, is an increase of 50 GeV for the three lightest masses 
and 60 GeV for the squark! If the threshold measurement is included, the effect is reduced 
to 30-40 GeV for the lighter masses, depending on the solution type, and 40-55 GeV for 
the squark. The increase is still surprisingly large and represents a serious threat to the 
applicability of the method in the case where the nominal masses happen to sit near a 
border. (Not to be confused with the previously mentioned border effect.) 
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It may however be possible to discover afterwards that such a mismeasurement has or 
may well have been made. If the masses obtained are used to generate the corresponding 
theory distributions, one may spot the problem. A quick test showed that in our case the 
m <j/(high) theory distribution as generated with the obtained masses, did have a small foot. 
This is then a sign that the original Tn'^^^s measurement may have been wrong, since no 
such foot was seen and included in the fit. This test probably has quite wide applicability. 
The reason is that even though, as we have seen, the masses may change a lot in the 
case of a mismeasured n^^g^y mass differences are not so far off the nominal values. 
The main change is therefore in the overall scale. For theory distributions, the shapes are 
independent of the overall scale. Only mass ratios are relevant. The theory shapes based 
on the masses obtained from the partly mismeasured endpoints, will therefore usually not 
deviate too much from the nominal shapes. 

If one is certain that a mismeasurement has been made, either from the test described 
above or from other considerations, not only does one know that the w^^gh) measurement 
is not to be trusted. There is also a way to pursue the matter further. It relies on the fact 
that the vertical fall we see and measure, is the endpoint of the m q i n distribution. We have 
the explicit expression for this endpoint, and so we can just replace the w^S gh ^ expression 
and redo the least square minimisation. 



At SPS la (f3) many new aspects of the endpoint method have emerged. For a large 
fraction of the experiments a total of three minima were competing for a given set of 
endpoint values, some returning masses near the nominal values, some returning much 
higher masses. A generic characteristic is the skewness of the ensemble distributions. The 
problems and the causes of the asymmetry have been discussed, but a way to handle the 
asymmetry has not been proposed. The ensemble distributions are much broader than 
at SPS la (a). Mass differences are however quite precisely determined. The problem of 
mismeasuring w^^gh)) which will be a danger in a large fraction of mass space, has been 
discussed and a solution proposed. Throughout, we have also seen that the inclusion of a 
fifth measurement in general improves on the results. 



7. Linear Collider inputs 



As compared with the LHC, the Linear Collider |4q , 49 1 will more directly provide very 
precise measurements of the lightest neutralino mass (and possibly also the second lightest). 
Thus, the LC input fixes the scale, whereas the cascade decays primarily provide mass 
differences. This kind of input, if it becomes available during the analysis of LHC data, 
will have a dramatic influence on the over-all analysis [27]. 

To each LHC experiment a corresponding Linear Collider experiment is considered. 
The Linear Collider measurement of the LSP, m~Jp, was for every experiment picked ran- 

Xi 

domly from a Gaussian distribution of mean equal to m-° m , and standard deviation set 

Xi 

by the expected uncertainty, o"^ c . We used <?m_ = 0.05 GeV, in accordance with [48]. 
Such a small uncertainty practically means fixing the mass at the nominal value. Still, for 
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completeness the measurement was appropriately included in the least square minimisation 



by adding the term [(m^o — ?™~ci ) / ^ml ] 2 to the E function, Eq. (| 
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Table 9: Number of minima for SPS la (/3) with LC input. 
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96.05 
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142.97 0.29 
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176.82 


176.82 0.17 




537.25 


537.2 2.5 
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491.92 


492.1 11.7 



Table 10: The SPS la (a) masses 
with LC input. All values in GeV. Re- 
gion (1,2) solutions now occur only in 
~ 1% of the cases and are left out. 



For (a) the fixing of m^o reduces the occurrences 
of multiple minima to the per mille level for any usable 
minima, AS < 3. In nearly all cases, 98-99%, it is 
the home region (1,1) minimum which survives. This 
seems reasonable since without the LC measurement 
the (1,2) minimum has (m^o) at some 10 GeV below 
the nominal value, see Table ^. 

At (/?) the high-mass minimum (1,1) is absent, 
disfavoured as it is by the Linear Collider measure- 
ment. The number of (1,3) minima has increased, 
which is reflected in an increase of the two-solution 
case, compare Tables fj] and |j| For small AS cuts 
this increase is not drastic. Fixing m^o thus does not 

X-l 

help us to uniquely determine one minimum, contrary to what one might have expected. 

The Tables [l0|-[Tl| show the mean and root-mean-square values of the ensemble dis- 
tributions at (a) and {(3). For both SUSY scenarios the ensemble means fall very close 
to the nominal masses, even for minima not situated in the nominal home region. The 
uncertainty on the scale, which for the LHC alone is the main contribution to the spread of 
the ensemble distributions, is set to zero by the LC measurement. The root-mean-square 
values are therefore strongly reduced. Without Linear Collider measurements, the mass 
differences were more accurately determined in that they were less dependent on the mass 
scale. A comparison of the root-mean-square values of the masses in Tables 0-|ll| with 
the root-mean-square values of the mass differences in Tables pHq shows to which extent 



mass differences are scale independent variables. When a Linear Collider measurement is 
available, the mass differences no longer out-perform the masses themselves in terms of 
precision. With the fixing of the scale also the skewness of the distributions has vanished. 



8. Conclusions 

We have investigated the measurement of supersymmetric masses from the decay chain 
Q ~^ X2Q ~^ IrIq ~^ Xv^l-i m the Snowmass mSUGRA scenario SPS la. Since the lightest 
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826.29 


826.1 6.3 


825.9 5.8 


828.6 5.5 



Table 11: The SPS la (/3) masses with LC input. Nominal, ensemble averages, (m), and root- 
mean-square deviations from the mean, a, are all in GeV. 

neutralino Xi is the LSP in most interesting mSUGRA scenarios, it will escape detection 
and only the quark and two leptons are available for the construction of invariant mass 
distributions. Nevertheless, the kinematic endpoints of these distributions have a well 
defined dependence on the masses of the particles in the decay chain and their measurement 
allows the extraction of the masses either by analytic inversion or numerical fit. The 
analytic expressions for the endpoints in terms of the masses were confirmed and presented 
together with their analytic inversions. 

In order to measure the endpoints of the invariant mass distributions pertaining to the 
chosen decay chain, one must have the correct mass hierarchies for the decay chain and a 
large enough cross-section. To ensure that this decay chain could be used over a wide range 
of scenarios, we performed a scan over the SUSY parameters. We found that as long as mo 
was not too large in comparison to mu2, a large proportion of the allowed parameter space 
would display the correct mass hierarchy. Furthermore, on examination of the sparticle 
production cross-sections and decay branching ratios we found that a large cross-section 
for the decay was available over much of this region. The Snowmass mSUGRA SPS la 
line/point falls into this region and is a good candidate for study. However, we noted that 
the cross-section for the decay chain is particularly high for the SPS la point, and it is 
instructive to examine a second point on the SPS la line with a less optimistic cross-section. 
We have denoted this new point as SPS la (f3) while the original point became SPS la (a). 

The LHC measurements of the endpoints were simulated using PYTHIA and ATL- 
FAST. Hard kinematic cuts remove practically all Standard Model backgrounds, except tt. 
Up to statistical fluctuations the powerful different-flavour subtraction then cancels the re- 
maining tt as well as any lepton-uncorrelated background from other SUSY channels. The 
resulting distributions are however contaminated by lepton-correlated background from 
x[>'s not taking part in the decay chain under study and combinatorial background from 
choosing the wrong jet. The inconsistency cut was shown to address the latter part with 
great efficiency, giving distributions much closer to the theoretical ones. Also mixed events 
were studied, revealing their potential to describe the background. More study is however 
needed. 

The endpoints were found by simplistic fitting of the edges, usually with a straight 
line together with a reasonable background estimate. A Gaussian smearing was sometimes 
included as a first approximation of the various smearing effects which take place. The 
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statistical precision of the edge was sought rather than an accurate determination of the 
endpoint. Still, the endpoints were seen to be in reasonably good agreement with the 
nominal values. However, the fitting procedure is clearly an area for improvement. On one 
hand, more realistic study of how the detector affects the distributions, and in particular 
the end regions, is called for. On the other hand, further study of the many realisations 
of the theory curves seems necessary. It is important to find good fit functions for the 
signal. A central part of such a programme is the incorporation of multiple squark masses 
at different rates. At a less ambitious scale the theory distributions should be studied 
for sheer acquaintance. The importance of such an awareness was demonstrated for the 
m qi(high) measurement at (/?). 

In order to turn the endpoint measurements into particle masses, and understand 
the resulting errors of this procedure, we considered an ensemble of 10,000 'gedanken' 
experiments. For each experiment a numerical fit for the particle masses was performed, 
using the method of least squares, thus appropriately handling the correlation between 
measurements due to the common jet energy scale error. Where available the analytic 
expressions for the masses in terms of the endpoints were used to provide starting points 
for the fits. 

The least squares function was found to often have two or even three minima of com- 
parable importance, a consequence of the multiple realisations for many of the endpoints. 
Without the threshold measurement there are for both scenarios usually two equally good 
minima, one in the correct region and one in another region, giving different masses. When 
the threshold endpoint is added, the minimum in the correct region is usually preferred. 
Still, in a noticeable fraction of experiments there will be more than one solution. Due to 
less precise measurements this applies more to {(5) than to (a). 

At (a) the minima of the correct region give masses very close to the nominal ones. The 
other (incorrect) region gives masses some 10-15 GeV lower. The ensemble distributions 
are symmetric. At ((3) there is one high-mass solution. The more precisely the threshold 
endpoint is determined, the less important this false minimum becomes. The low-mass 
solutions, one or two, are closer to the nominal values, but the distributions are skewed. 
This is a combined effect of the large endpoint uncertainties and the so-called border effect. 

The obtained masses of the three lightest particles are found to be very strongly related. 
Furthermore it was seen that mass differences are better variables, in the sense that they 
are less correlated than the masses themselves. Due to the form of the endpoint expressions, 
the LHC will measure mass differences at high precision, but leave the overall scale less 
certain. A Linear Collider measurement of the LSP mass effectively sets the scale, which is 
why the precision of the masses improve drastically when the LHC and the Linear Collider 
measurements are combined. 
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